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DIFFEOMORPHIC APPROXIMATION OF 
ipi 1 PLANAR SOBOLEV HOMEOMORPHISMS 


STANISLAV HENCL AND ALDO PRATELLI 

Abstract. Let H C be a domain and let / G be a homeomorphism 

(between fi and /(n)). Then there exists a sequence of smooth diffeomorphisms fk 
converging to / in and uniformly. 


1. Introduction 

The general problem of finding snitable approximations of homeomorphisms / : M"' D 
n — f{fl) C M" with piecewise affine homeomorphisms has a long history. As far as 
we know, in the simplest non-trivial setting (i.e. n = 2, approximations in the L°°-norm) 
the problem was solved by Rado PI. Dne to its fnndamental importance in geometric 
topology, the problem of hnding piecewise affine homeomorphic approximations in the L°°- 
norm and dimensions n > 2 was deeply investigated in the ’50s and ’60s. In particnlar, it 
was solved by Moise [21] and Bing [7] in the case n = 3 (see also the snrvey book [25]). 
while for contractible spaces of dimension n > 5 the resnlt follows from theorems of 
Connell [2], Bing [H], Kirby [20] and Kirby, Siebenmann and Wall [21] (for a proof see, 
e.g.. Rushing [3T] or Luukkainen [22] )• Finally, twenty years later, while studying the 
class of quasi-conformal varietes, Donaldson and Sullivan [T2] proved that the result is 
false in dimension 4 in the context of manifolds. 

Once completely solved in the uniform sense, the approximation problem suddenly 
became of interest again in a completely different context, namely, for variational models 
in nonlinear elasticity. Let us briefly explain why. In the setting of nonlinear elasticity (see 
for instance the pioneering work by Ball [3]), one is led to study existence and regularity 
properties of minimizers of energy functionals of the form 

/(/)= / W{Df)dx, (1.1) 

Jn 

where / : M” D D —)■ A C M"' (n = 2, 3) models the deformation of a homogeneous elastic 
material with respect to a reference conhguration D and prescribed boundary values, while 
W : —)■ M is the stored-energy functional. In order for the model to be physically 

relevant, as pointed out by Ball in [H |5], one has to require that u is a homeomorphism 
-this corresponds to the non-impenetrability of the material- and that 

W{A)^+oo asdetA—)-0, IF(A) =-|-oo if det A < 0 . (1.2) 
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The first condition in (1.2) prevents from too high compressions of the elastic body, while 
the latter guarantees that the orientation is preserved. 

Another property of W that appears naturally in many problems of nonlinear elasticity 
is the quasiconvexity (see for instance BD- Unfortunately, no general existence result is 
known under condition (1.2), not even if the quasiconvexity assumption is added: one has 
either to drop condition (1.2) and impose p-growth conditions on W (see m^) , or to 
require that W is polyconvex and that some coercivity conditions are satished (see mm)- 
Moreover, also in the cases in which the existence of minimizers is known, very little 
is known about their regularity. 

As pointed out by Ball in [U |5] (who ascribes the question to Evans 1131). an important 
issue toward the understanding of the regularity of the minimizers in this setting (i.e., W 
quasiconvex and satisfying (1.2)) would be to show the existence of minimizing sequences 
given by piecewise affine homeomorphisms or by diffemorphisms. In particular, a hrst 
step would be to prove that any homeomorphism u G M""), p G [1, +C)o), can be 

approximated in by piecewise affine ones or smooth ones. One of the main reasons 
why one should want to do that, is that the usual approach for proving regularity is to test 
the weak equation or the variation formulation by the solution itself; but unfortunately, 
in general this makes no sense unless some apriori regularity of the solution is known. 
Therefore, it would be convenient to test the equation with a smooth test mapping in the 
given class which is close to the given homeomorphism. More in general, a result saying 
that one can approximate a given homeomorphism with a sequence of smooth (or piecewise 
affine) homeomorphisms would be extremely useful, because it would signihcantly simplify 
many other known proofs, and it would easily lead to stronger new results. It is important 
to mention here that the choice of dimension n = 2, 3 is not only motivated by the physical 
model, but also by the fact that the approximation is false in dimension n > 4, as shown 
in the very recent paper mi. 

However, the Ending of difieomorphisms near a given homeomorphism is not an easy 
task, as the usual approximation techniques like mollification or Lipschitz extension using 
maximal operator destroy, in general, the injectivity. And on the other hand, we need of 
course to approximate our homeomorphism not just with smooth maps, but with smooth 
homeomorphisms (otherwise the approximating sequence would be not even admissible 
for the original problem). 

Few words have to be said about the choice of the required property for the approxi¬ 
mating sequence, namely, either smooth or piecewise affine. Actually, both results would 
be interesting in different contexts. Luckily, the two things are equivalent: more precisely, 
it is clear that an approximation with difieomorphisms easily generates another approxi¬ 
mation with piecewise affine homeomorphisms; the converse is not immediate but, at least 
in the plane, it is anyway known (see m- Therefore, one can approximate in either of 
the two ways, and the other one automatically follows (for instance, in this paper we will 
look only for piecewise affine approximations). It is important to clarify a point: whenever 
we say that a map is piecewise affine, we mean that there is a locally finite triangulation 
of such that the map is affine on every triangle. It is actually possible to find finite 
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triangulations whenever this makes sense; but for instance, if hi is not a polygon, then 
the triangles must obviously become smaller and smaller near the boundary, so a finite 
triangulation does clearly not exist. 

Let us now describe the results which are known in the literature in this direction. The 
first ones were obtained in 2009 by Mora-Corral |2S] (for planar bi-Lipschitz mappings that 
are smooth outside a finite set) and by Bellido and Mora-Corral [B], in which they prove 
that if u, u~^ G (7°’“ for some a G (0,1], then one can find piecewise affine approximations 
n of M in (7°’^, where {3 G (0, a) depends only on a. 

More recently, Iwaniec, Kovalev and Onninen [18] almost completely solved the ap¬ 
proximation problem of planar Sobolev homeomorphisms, proving that any homeomor- 
phism / G M^), for any 1 < p < -|-cx), can be approximated by diffeomorphisms fe 

in the norm (improving the previous result, of the same authors, for the case p = 2, 
see jini). 

Later on, it was shown by Daneri and Pratelli in [101 and [H] that any planar bi- 
Lipschitz mapping / can be approximated by diffeomorphisms fk such that fk converge 
to / in norm and simultaneously converge to f~^ in W^’^, giving the first result 
in which also the distance of the inverse mappings is approximated. 

The goal of the present paper is to prove the approximation of planar homeo- 
morphism in the sense, so basically dealing with the important case p = 1 which 
was left out in [TH|. In particular, our main result is the following. 

Theorem 1.1. Let 12 C &e an open set and f G hL^’^(12, M^) be a homeomorphism. For 
every e > 0 there is a smooth diffeomorphism (as well as a countably -but locally finitely- 
piecewise affine homeomorphism) fs G 1K^’^(12,M^) such that \\fe—f\\w^v +II/.-/II [^oo < e. 
If in addition f is continuous up to the boundary o/12, then the same holds for fe, and 
f^ = f on dQ. 


Actually, our piecewise affine functions fe will be globally finitely piecewise affine, 
thus also bi-Lipschitz, as soon as 12 is a polygon and / is piecewise affine on dfl, see 
Theorem 4.20 If even just one of these two conditions does not hold, then this is clearly 


impossible (see Remark 4.21). 


We conclude the introduction with a short comparison between the techniques of this 
paper and those of the other papers discussed above. The proofs in |26l |6| are based on a 
clever refinement of the supremum norm approximation of Moise |23|, while the approach 
of [H] and of the other contributions of the same authors makes use of the identification 
~ C and involves coordinate-wise p-harmonic functions. The techniques of the present 
paper are completely different with respect to them; basically, our proof is constructive, 
and it is based on an explicit subdivision of the domain / which depends on the Lebesgue 
points of Df. 

Our techniques resembles the basic ideas of [10] and m, and we will also use some 
of the tools introduced there, but there are some extremely important differences. More 
precisely, on one hand in ra and HH one had to approximate at the same time Df and 
Df~^ while here we need only to approximate Df, and this is a deep simplification. But 
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on the other hand, in this paper we look for a sharp estimate, that is, / is only in 
and we want an approximation exactly in thus we have not so regular maps and 

we cannot lose sharpness of the power anywhere, while in [10] and m the maps were 
much better, namely bi-Lipschitz, and in several steps the sharpness of the power was lost. 
Roughly speaking, we can say that the most difficult steps of [10] and m correspond to 
much simpler steps here, and vice versa. 

1.1. Brief description of the proof. In this section we outline the basic plan of our 
proof, to underline the main steps and help the reading of the construction. We remind 
the reader that our aim is to find an approximation done with piecewise affine homeomor- 
phisms, and then the existence of an approximation with smooth diffeomorphisms will 
eventually immediately follow applying the result of iza. 

First of all, we will divide our domain into some locally finite grid of small squares, 
these squares becoming maybe smaller and smaller close to dfl. We will then consider 
separately the “good” squares, and the “bad” ones. More precisely, a square S (c, r) in the 
grid will be called good if / can be well approximated by a linear mapping f{c)+M{x — c) 
there, where M coincides with Df in some Lebesgue point close to c: in particular, we 
will need that j-s{c,r) \^f ~ is small enough. Since almost every point of hi is a 

Lebesgue point for Df, we will be able to deduce that, up to consider a sufficiently hne 
grid, the area covered by the good squares is as close as we wish to the total area of hi. 

Moreover, up to a slight modihcation of the value of / on the boundary of the squares, 
we will reduce ourselves to the case that 

\Df\<K^\Df\, (1.3) 

JdS Js 

where K is a big, but fixed, constant. 

We will then define an approximation of / (which will eventually become f^) on the 
grid: on the boundary of each square, we will hnd a piecewise linear approximation of 
/, very close to /, in such a way that these approximations on the whole grid remain 
one-to-one (we do not have just to take care of the approximation on a single square, but 
also check that the different approximations coincide on the common sides, and that they 
do not overlap with each other). Of course this will be much easier for the good squares, 
since on a whole good square / is already almost affine, and more complicated for the 
bad squares. We will do our approximation g in such a way that, for any bad square S, 

\Dg\<K f \Df\ 

JdS 

The next step is then to extend the piecewise linear maps to the interior of each square; a 
good thing is that, being g already defined on the grid in a one-to-one way, the extension 
inside each square is completely independent with what happens on the other squares. 
The rough idea to do so is that on good squares we can obtain very good estimates, while 
in bad squares we can get only bad estimates; but since the total area of the bad squares 
is arbitrarily small, in the end everything will work. 

The first tool which we will need, presented in Section says that any piecewise 
linear map g defined on the boundary of a square S can be extended to a piecewise affine 


(1.4) 
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homeomorphism h in the interior of S in such a way that 

-f-\Dh\<K-f- \Dg\. (1.5) 

JS JdS 

This construction is done hrst by choosing many points on the boundary of the square; 
then, for any two of these points, say x and y, we will select the shortest path joining 
g{x) with g{y) remaining inside the portion of having as boundary. Using in a 

careful way these shortest paths, we will then eventually obtain the definition of h such 


that (1.5) holds true. This estimate, together with (1.3) and (1.4), readily implies that 


for every bad square S one has 


\Dh\<K ( \Df\- 
Js 


, , , . . , ( 1 - 6 ) 
Is Js 

we have then just to take a very hne grid, so that a very small portion his of hi is covered 
with bad squares, and hence we will have 

/ \Df-Dh\<j \Df\ + \Dh\<{K + l) j \Df\<e. 

It remains then to consider the good squares, and here we will have to be extremely precise. 
As already said, around every good square S the map / is very close to being affine, hence 
the image of S is very close to a parallelogram; therefore, there is no problem unless this 
parallelogram degenerates. Let us be more precise: for all the squares corresponding 
to a matrix M with strictly positive determinant (hence, the parallelogram does not 
degenerate), the extension inside S is trivial; it is enough to divide the square in two 
triangles and consider on each triangle the affine map which equals / on the three vertices. 
By construction, we will easily see that this works perfectly. 

The good squares corresponding to M = 0 are a problem only in principle: indeed. 


we can treat them as bad squares. The estimate (1.6) says that this gives a cost of a big 


constant K times the total integral of Df on those squares; however, since they are good 
squares and the corresponding matrix is M = 0, by definition the integral of Df will be 
extremely small, and everything will work. 

The hard problem, instead, is for good squares for which M 7 ^ 0 but det M = 0: these 
correspond to degenerate parallelograms, and we have to treat them carefully because 
these squares can cover a large portion of hi: in fact, recall that the set {x : \Df{x)\ 7 ^ 
0 and Jf{x) = 0} can have positive or even full measure for a Sobolev homeomorphism 
(see USD- Section is devoted to build the extension for this specific case, which is 
somehow similar to the one with the shortest paths described above. The big difference 
with that case is, on one hand, that this time we are in a good square, hence very close to 
a Lebesgue point for Df, and this helps even for the degenerate case. But on the other 


hand, this time we are not satisfied with an estimate like (1.5), where a big constant K 


appears, and we need instead an approximation h which is very close to the original /. 
This extension procedure will be the most delicate step of the construction. 

The construction of the proof, divided in its several steps, is done in the last Section]^ 
Basically, putting together all the ingredients described above, the proof will then be 
concluded for what concerns the existence of the piecewise affine approximation; the 
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existence of the smooth approximation will then follow thanks to the result of [ 2 Z] , while 
the claim about the boundary values will be easily deduced by the whole construction. 


1.2. Preliminaries and notation. In this section we shortly list the basic notation that 
will be used throughout the paper. By S{c,r) we denote the square centered at c, side 
length 2 r and sides parallel to the coordinate axis, while Sq = {(x, y) G : |x| + \y\ < 1 } 
is the “rotated square”, that we use only in Section]^ Similarly, B{c, r) is the ball centered 
at c with radius r. 

The points in the domain hi will be usually denoted by capital letters, such as A, B 
and so on, while points in the image /(hi) will be always denoted by bold capital letters, 
such as A, B and similar. To shorten the notation and help the reader, whenever we use 
the same letter A for a point in the domain and A (in bold) for a point in the target, this 
always means that A is the image of A under the mapping that we are considering in that 
moment. By AB (resp., AB) we denote the segment between the points A and B (resp. 
A and B). The length of this segment is denoted as 'H}{AB), or as AB^ while 
is the length of a curve 7 . With the notation AB (or AB) we will denote a particular 
path between A and B (or A and B), whose length will be then T-i^iAB), or 'H^(AS); 
we will use this notation only when it is clear what is the path we are referring to (often 
this will be a shortest path between the points). Given three non-aligned points A, B, C 
(or A, B, C), we will denote by ABC (or ABC) the angle in (0, vr) between them, and 
by ABC (or ABC) the triangle having them as vertices. 

We will denote the (modulus of the) horizontal and vertical derivatives of any mapping 

/ = (/i,/2):M2^M2 ^3 


Dif\ = 



^ 2/1 




Analogously, the derivatives of the components fi and /2 are written as 




dh 

dx 


D2fl 


dh 

dy ’ 


Dlf2 


df2 

dx 


D2h 


df2 

dy ' 


Whenever a continuous function g is dehned on some curve 7 (usually, 7 will simply be 
the boundary of a square) we will denote by r(f) the tangent vector to 7 at f G 7 , and 
by Dg{t) the derivative of 5 ^ at f in the direction of T{t). With a small abuse of notation, 
even if the derivative is not necessarily defined, we will write \Dg{t) \ (BH}{t) to denote 
the length of the curve gh)'. notice that the latter length is always well-defined, possibly 
-|-oo, and it actually coincides with hs \Dg{t)\ as soon as this is defined. Finally, notice 
that, if a function / is affine on a square 5, being Df = M for some matrix M, and we 
call g the restriction of / to dS, then Dg(t) = M ■ rit) for any t G dS. 

The letter K will always be used to denote a large purely geometrical constant, not 
depending on anything; we will not modify the letter, even if the constant may always 
increase from line to line. For the sake of simplicity (and since the precise value of K 
does not play any role) we do not explicitely calculate the value of this constant. 
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2. Extension from the boundary of the square 


This section is entirely devoted to show the result below about the extension of a map 
from the boundary of the square to the whole interior. 

Theorem 2.1. Let g : OSq —)■ be a piecewise linear and one-to-one function. There is 

a finitely piecewise affine homeomorphism h \ Sq ^ such that h = g on dSo, and 


/-So 


\Dh{x)\dx<K / \Dg(t)\d'H^(t). 


'dSo 


( 2 . 1 ) 


Proof. The construction of the map h is quite long and technical, and hence we subdivide 
it in several steps. 

Step 1. Choice of good corners, so that (2.2) holds. 

For our construction, we will need to assume that \Dg\ does not concentrate too much 
around the corners; more precisely, we will need that 

[ \Dg\dH^<Kr[ \Dg\dH^ for all r G (0,1), i e {1, 2} , (2.2) 

JB{Vi,r)r\dSo JdSa 


being Vi = (0, —1) and V 2 = (0,1). It is quite easy to achieve that: in fact, it is enough 
to hnd two opposite points Pi, P 2 & OSq such that 


[ \Dg\ <6r[ \Dg\ for all r G (0, y/2), i G {1, 2} , 

J!3(Pi,r)ndSo JdSo 


(2.3) 


because then we can apply a bi-Lipschitz transformation (with bi-Lipschitz constant in¬ 
dependent from Pi and P 2 ) which moves the points Pi and P 2 on the vertices Vi and V 2 , 
and get (2.2). And in turn, to obtain (2.3), we notice that every point of dSo is a possible 
choice for Pi or P 2 unless it is, or its opposite point is, in the set 


^ := P G : 3r G (0,1) : 


JB{P,r)ndSo 


\Dg\ >6r [ \Dg\ 
JdSo 


By a Vitali covering argument, we can cover A with countably many balls B{Pi, Srf such 
that every Pi is in A, and the corresponding sets P(Pj, rf ft are as in the dehnition 
of A and are pairwise disjoint. Therefore, we can calculate 

nHA) <Y.6r, < y: 'f*"' < 1 , 

i i Jaso 

and since PL^{dSo) = 4^/2 it clearly follows that two opposite points both in dSo \.4, exist 
and then satisfy (2.3), as required. 

Step 2. Definition of the grid on dSo, and of the paths 7 *. 

To dehne our map h, we will make use of a hne grid made by horizontal segments 
in So. More precisely, we will take several (but hnitely many) distinct points = 

= (0,1) in all with non-positive hrst coordinate Aj and 
with second coordinate A^ increasing, with respect to i, from —1 to 1; on the opposite 
side, we will take the corresponding points P* = {—AfAf), so that the segments A*P* 
are horizontal. 

The way to choose our points is simple: since g is piecewise linear, we can take the 
points in such a way that g is linear on every segment A^A^^^, as well as in every 
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Since this property is of course not destroyed if we add more points Ad (as long as we also 
add the corresponding points 5*, of course), we are allowed to add more points during 
the construction, of course taking care to add only hnitely many: we will do this a hrst 
time in few lines, and then also later. 

From now on, we will call «So the bounded component of \ g{dSQ), which is a 
polygon because g is piecewise linear; notice that the map h that we want to construct 
must be a homeomorphism between So and «So. Then, for any 0 < i < k, we dehne 7 * 
the shortest path which connects A* and inside the closure of «So (this shortest path 
is unique, as we will show in Step 3). Notice that, since «So is a polygon, every 7 * is 
piecewise linear, and any junction between two consecutive linear pieces is in dSo. 

Up to add one more point between and A^ (plus the corresponding one on the 
right part), we can suppose that 7 ^ is either a segment between and B^, and this 
happens if and only if the angle between A° and B^ which goes inside «So is smaller 
than TT, or it is done by the union of the two segments A^A° and A^B^, thus it entirely 
lies on OSq. We do the same between A^~^ and A^. 

Step 3. Uniqueness of the shortest paths. 

Let Q be a simply connected closed planar domain with polygonal boundary. We briefly 
recall the proof of the well-known fact that, for any two points in Q -not necessarily on 
the boundary- there is a unique shortest path inside Q. Since the existence is obvious, 
we just have to check the uniqueness. 

If the claim were not true, there would be two points A, B E Q and two shortest paths 
Ti and T 2 between A and B inside Q, such that ti and T 2 meet only at A and B. The 
union of the two paths is then a polygon, say with n sides. The sum of the internal angles 
of this polygon is 7r(n — 2), and thus there must be a vertex of the polygon, different from 
A and B and thus inside one of the shortest paths, corresponding to an angle strictly less 
than TT. Since the interior of the polygon is entirely in the interior of Q, this is of course 
impossible, because cutting around that vertex would strictly shorten the length of the 
path, against the minimality. 

Step 4- The path is above 7 *, and definition of jI, 

For two curves 7 and 7 inside «So and with endpoints in OSq, we say that “7 is above 
7 ” if 7 does not intersect the interior of the (possibly disconnected) subset of «So whose 
boundary is the union between 7 and the path on dSo connecting the endpoints of 7 and 
containing Afi = g{A^). We want to show that, for any 0 < i < /c — 1, the path 7 *+^ is 
above 7 *. 

To show that, assume that two points P and Q belong to both the paths 7 * and 7 *+L 
Then, the two restrictions of 7 * and from P to Q are two shortest paths, and by Step 3 
we derive that 7 * and coincide between P and Q. As an immediate consequence of 
this observation, we get that is above 7 * as claimed. 

Another immediate consequence is the following: the intersection between 7 * and 
is always a connected subpath, possibly empty. If it is not empty, and then it is a path 
PQ, we will subdivide both 7 * and in three parts, writing 

7 ' = yj U 7 * U 7 * , 7 '+^ = yj+i U 7 *+i U 7 *+^, 
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where 7 ^ (resp. 7 ^"'"^) is the hrst part, from A* to P (resp., from to P); 72 (resp. 
72 '''^) is the second part, from P to Q (thus the common part, and 72 = 72 '*'^); and 73 
(resp., 73 ’'"^) is the third and last part, from Q to P* (resp., from Q to If 7 * and 

7 *'''^ have empty intersection, then we simply set 7 ^ = 7 * and 71 ^^ = 7 *^^, letting 72 , 73 , 
72 '*'^ and 73 ^^ be empty paths. The situation is depicted in Figure where the common 
part 72 = 72 ’'"^ is done by two segments, one on dSo and the other one in the interior of 
5o. 



Figure 1. The paths 7 * and 7 *’''^ in Step 4. 

Notice that this subdivision of a path does not depend only on the path itself, but 
also on the other path that we are considering; in other words, the subdivision of the path 
7 -^ done when i = j, and then considering the possible common part between 7 -^ and 7 -^’'"^, 
does not need to coincide with the subdivision of the same path done when i = j — 1, and 
then considering the possible common part between 7 -^ and 7 -^“^. 

Step 5. Convexity of the polygon having boundary U A*’''^P. 

Let us call P the last point of the path 71 ’''^; hence P is the hrst common point with 
7 *, if 7 * and 7 *’''^ have a non-empty intersection, while otherwise P = We claim 

that the polygon having U A*'''^P as its boundary is convex (notice that in principle 
the curve 7 ^^^ and the segment A*’'"^P could have other intersection points in addition 
to A*^^ and P). We start assuming that 7 ^ 0, at the end of this step we will then 
consider the other case. 

If yj"*"^ is a single point, or just a segment, then the claim is emptily true, and the 
convex polygon is degenerate. Let us assume then that yj^^ is done at least by two 
affine pieces, and assume also, just to hx the ideas, that the direction of the oriented 
segment A^A*"*"^ is 7 r/ 2 , as in Figure left. Call then P So the polygon having, 
as boundary, the Jordan curve yj U yj’'"^ U A*^^AL The same argument as in Step 3 
immediately ensures that, for any vertex of (i.e., any junction point between two 
consecutive linear pieces of yj"'"^), the angle pointing inside <So (hence in particular inside 
P) is bigger than tt. By construction, and recalling Step 4, we get also that none of 
these points can belong to the curve in dSo connecting A® and P® and containing A°, 
since such a point should necessarily belong also to 7 ®, against the dehnition of P. Of 
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course, this already “suggests” that our convexity claim is true, but observe that the proof 
is still not over, since in principle could be some spiral-like curve connecting 
with P. To conclude the proof, for any vertex of (except P) consider the range of 
directions pointing toward the interior of T>\ for instance, the range associated to in 
the situation of Figure left, is done by the angles between — 7 r /2 and —vr/S. We claim 
that, for each vertex of the curve 71 ’'"^, this range cannot contain the angle -|-7r/2: observe 
that this will immediately imply the required convexity. 

Assume then by contradiction that this claim is false, and let Q be the hrst vertex 
of 7 ];’''^ having 7r/2 in its range of directions; by a trivial perturbation argument we can 
assume that 7r/2 is in the interior of this range, and then the vertical line passing through 
Q is in the interior of T> for a while, both above and below Q itself. Call then, as in 
Figure left, Q~ and the hrst points of this line, respectively below and above Q, 
which are on dP. Since the segment Q Q^ is parallel to A^A*"*"^, each of these points 
must belong either to 7 ^ or to 71 '''^. Observe now that, if 7 is a shortest path in Sq between 
its extremes, it is also a shortest path in Sq between any pair of its points. In particular, 
if the line connecting two points of 7 is entirely in the closure of Sq, then 7 must be the 
segment between these two points. This immediately imply that none of the points Q 
and can belong to 71 '''^, because otherwise should be a segment between that 
point and Q; as a consequence, both the points Q~ and must belong to 7 ^, but this is 
also impossible because then yj should be the segment between them. The contradiction 
shows the claim, and then we have obtained the required convexity. Of course, the very 
same argument works for the polygon having boundary yg"*"^ U being this time 

P the hrst point of 73 ^^, and everything also works for the polygons around the path 7 ® 
instead of 7 ®+^. 


^i+i 



Q 


Q 



Figure 2. Construction in Step 5. 

Let us then consider the case when yg"*"^ = 0, that is, the case when 7 ® and y®"*"^ 
are disjoint: this situation is depicted in Figure right. This time, we dehne P the 
polygon having, as boundary, the Jordan curve 7 ® U .B®.B®’'"^ 7 ®’''^ U A®^^ A®. The very same 
argument as in the hrst case ensures again that every vertex of y®"*"^ has angle bigger than 
71 in the direction inside X>; as a consequence, the required convexity follows, as before, if 
the range of every vertex of y®"*"^ does not contain the direction -|-7r/2. 
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However, this time it is not impossible that a vertex Q of has 7 r /2 in its range. 
Let then, as before, Q be the hrst vertex (if any) with this property, and let G be 
as before. As already noticed, none of the points can be in and at most one in 
7 *. Hence, the only possibility is that one point is in 7 *, and the other one in A 

simple topological argument ensures that Q~ must be in 7 * and in B'^B'^^^. Indeed, 
consider the path, contained in dT> and not containing which connects Q and 

Q~] together with the segment Q Q, this is a Jordan curve, and then it can not intersect 
the other path in which contains in particular, it must contain Q^, and 

it readily follows, as claimed, that Q G 7 *, G BiB^^^. The very same topological 
argument ensures also that is the “left” vertex (that is, the one in the direction 

A*A*^i) of the segment B^B^~^^, and B^ is the “right” one, as in Figure]^ right. 

Let us now restrict our attention to the subset X>o of made by the polygon whose 
boundary is the part of 7 *+^ connecting Q to B^^^, plus the two segments and 

Q+Q. The same argument of the first half of this step ensures that the range of directions, 
toward the interior of X>o, corresponding to any vertex of 7 *+^ in can never contain 

the direction of the segment B^B^~^^, since otherwise a segment parallel to B^B^~^^ and 
contained in X>o should have both the endpoints in the segment QQ^, which is impossible. 

Finally, it is immediate to notice that this property of the directions, analogously as 
before, is enough to ensure the required convexity of the polygon having U as 

boundary. 

Step 6. Definition of the “vertical segments” and their length. 

In this step, we associate to any vertex P of the curve 7 *+^ a point (or many points) Q 
of the curve 7 ®, and vice versa. Every such segment PQ, which we will call “vertical”, 
will be contained in the closure of the polygon X> C <So having boundary A^A®"*"^ Uy®’*'^ U 
Uy®, any two vertical segments will have empty intersection, except possibly at a 
common endpoint, and the following estimate for the length of the vertical segments will 
hold, 

n\PQ) < max {h^(A'A'+^), H^(M®+^)} . (2.4) 


Let us give our definition distinguishing the possible cases, as in Step 5. 

First of all, consider the situation, depicted in Figure left, when 7® and 7®+^ have 
a non-empty intersection. In the common part 7® fl 7®+^ = 1 2 — 12^'^i ’wiH associate 

to any vertex P of y®"*"^ the same point Q = P, which is also in 7® by dehnition. The 
segment PQ is just a point, which is of course in the closure of T>, and the length is 
0, so that (2.4) of course holds. In the “left” part of the paths, instead, we will give 
the following simple dehnition. To any vertex P G yj"'"^, we associate the point Q G yj 
so that the segment PQ is parallel to A®’'"^A®: the existence and uniqueness of such a 
point, the validity of (2.4), and the fact that PQ is contained in the closure of T>, all 
come immediately from the convexity obtained in Step 5. We do the very same thing for 
the vertices of yj, and we argue completely similarly for the “right” part of the paths, 
of course taking segments parallel to B^^^B^, instead of A®’'"^A®. Then we have already 
completed our dehnition of the vertical segments, and the fact that any two such segments 
do not intersect is obvious from the construction. 
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_g*+i 




Figure 3. Construction in Step 6: the polygon «So (resp., T>) is light (resp., 
dark) coloured, and the “vertical segments” are dotted. 


Consider now the situation when 7® n 7®’*'^ = 0 , see Figure right. Without loss of 
generality we can think that, as in the Figure, the direction of is vertical, while 

the segment goes “toward left”. Let then S G and T G 7® be the two closest 

points such that the segment TS is vertical; notice that it is possible that S = A®"*"^ 
or that S = .B®"*"^, this makes no difference in our proof, even if the picture shows an 
example where S is in the interior of the curve y®"*"^. 

Let us now consider the subset T>o C T> whose boundary is given by the two segments 
A®A®^^ and ST, together with the parts of 7® (resp., 7®^^), connecting A® and T (resp., 
A®"*"^ and S). Again by the convexity result of Step 5, it is clear that at any point of y®"*"^ 
between A®"*"^ and S starts a vertical segment, whose interior is entirely contained in X>o, 
which ends in a point of 7® between A® and T. We define then in the obvious way the 
“vertical segments” inside 'Dq, which are in fact vertical. The validity of (2.4) is as usual 
obvious from the convexity. 

Consider now the half-line starting at S and parallel to B®^^B®. The choice of the 
points S and T, together with the convexity proved in Step 5, ensure that this half-line 
remains inside for a while, after the point S] therefore, the intersection of this half-line 
with X> is a segment ST^, and the point T’*' is on 7® by construction. Observe that 
coincides with T in the particular case when B®’'"^B® is parallel to A®’'"^A®, but otherwise 
it stays, as in the hgure, outside of X>o. The construction implies that all the half-lines, 
parallel to B®^^B® and starting at a point of y®"*"^ after S, remain in T> for a while and 
then intersect 7® at some point after T^. We use this observation to associate to any 
vertex of 7®^^ after S a point of 7® after T’*', and we call then “vertical segments” all the 
corresponding segments, which are actually not vertical but parallel to B^~^^B^. Finally, 
to every vertex of 7® between T and T^, if any, we associate always the point S. The 
validity of (2.4) for all the vertical segments is then again clear by the construction and 
by Step 5, and any two vertical segments have always empty intersection, unless in the 
case when they meet at S. This concludes the step. 
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From now on, we will always consider as “vertices” the points S, T and T^, even 
if they were not vertices in the sense of the piecewise linear curves. Moreover, for every 
vertex of 7 *, or of we will consider as “vertex” also the corresponding point in the 
other curve, which again could be or not be a vertex in the classical sense. Notice that 
in this way we are adding a finite number of new vertices and, as already pointed out 
before, it is always admissible to regard as “vertices” also finitely many new points in our 
curves. Summarizing, on the piecewise linear curve 7 ® we are considering as “vertices” 
all the actual vertices, plus some other new points. However, these “new points” have 
been selected working on the region between 7 ® and y®"*"^, and then they do not need to 
coincide with the “new points” selected by working on the region between 7 ®“^ and 7 ®. 
Step 7. Definition of h on Sq. 

We are now ready to define a function on So which extends g; for simplicity, we start now 
with the definition of a “temptative” function h, without taking care of the injectivity. 
The definitive function h will be obtained later. 

Recall that we have selected several horizontal segments A^B^, 1 < i < fc — 1, in the 
square Sq] the square is then divided in k — 2 “horizontal strips”, i.e. the regions between 
two consecutive horizontal segments, plus two triangles, the “top one” and 

the “bottom one” A^A^B^. 

We start defining the function h on the “1-dimensional skeleton”, that is, the union of 
dSo with all the horizontal segments A^B^: more precisely, we set h = g on the boundary 
dSo, while for every 1 < i < k — 1 we define h on the horizontal segment H®i?* as the 
piecewise linear function, parametrized at constant speed, whose image is the path 7 ®. 
Notice that, with this dehnition, h is continuous on the 1-skeleton. 

To extend h to the whole 5o, we can then argue separately on each of the horizontal 
strips of So, as well as on the top and bottom triangle. First, let us consider the bottom 
triangle A^A^B^\ thanks to the construction of Step 2 , we know that the path 7 ^ is 
either the segment A^B^, or the union of the two segments and A^B^. In the hrst 

case, we dehne h on the bottom triangle as the affine function extending the values on the 
boundary; in the second case, let P be the point of the segment A^B^ such that h{P) = A^, 
let us extend h as constantly A^ on the segment PA^, and let h be the (degenerate) affine 
function extending the values on the boundary on each of the two triangles A^PA^ and 
/PPB^. In the top triangle, we give of course the very same dehnition of h. 

Let us now consider the horizontal strip Di between H®i?® and and let us 

call the bounded region in <So having as boundary the closed curve y®"*"^ U B^^^B^ U 
7 ® U A''A''^^. In Step 6 , we have selected a hnite number of points on 7 ® and on y®"*"^, 
and we have called “vertical segments” the corresponding segments. More precisely, let 
us denote the points in H®+^R®+^ as Pq = H®+\ Pi, ..., Pm-i, Pm = P®"*"^, and the points 
in H®P® as Qo = Af Qi, . .., Qm-i, Qm = P*; as always, let us write Pj = h{Pj), and 
Qj = h{Qj). Keep in mind that each segment PjQj, whose interior is entirely contained 
in T>i, has been called a “vertical segment”, and notice that the points Pj and Qj are 
not necessarily all different: for instance, the point S of Figure]^ right, is the point Pj 
for three consecutive indices 0 < j < M. 
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We are finally in position to give the definition of h on the interior of each strip Vi (and, 
since h has been already dehned in the 1-skeleton and on the top and bottom triangle, 
this will conclude the present step). The strip Vi is the essentially disjoint union of the 
triangles PjPj+iQj and Pj+iQjQj+i for all 0 < j < M, and Vi is the essentially disjoint 
union of the corresponding triangles PjPj+iQj and Pj+iQjQ where the triangles 
in Vi (but not those in Vi) can be degenerate, in particular they are degenerate for the 
points in We dehne then h on Vi as the function which is affine on each of the 

above-mentioned triangles. Notice that, by construction, h is linear on each side PjPj+i 
and QjQj^i, hence this dehnition on Vi is a continuous extension of the dehnition on the 
1 -skeleton. 

Step 8. Estimate for \Dh\. 

In this and in the following step, we aim to estimate the integral of \Dh\ on iSq; in 
particular, in this step we will consider the bottom triangle (by symmetry, we 

will get an estimate valid also for the top triangle, of course), while in the next step we 
will consider the situation of the horizontal strips Vi. The aim of this step is to show the 
validity of the bound 

[ \Dh\<K [ \Dg\dH\ (2.5) 

Jaoa^B^ JdSo 

where as usual K denotes a purely geometric constant. By simplicity, let us call 


f :=n\A^A^) = n\A^B^). 


Recall that, on the bottom triangle, the function h has been dehned as an affine function, 
if the angle pointing inside «So, is smaller than tt -or, equivalently, if the curve 

7 ^ coincides with the segment A^B^- and as two degenerate affine functions on the two 
triangles A^PA^ and A^PB^ (being P as in Step 7) otherwise. Let us then estimate the 
norm of Dh on the bottom triangle in both cases. 

First of all, consider the non-degenerate case when h is a single affine function on 
the bottom triangle. In particular, the image of the segment A^A^ is the segment A^A^, 
while the image of the segment A^B^ is the segment A^B^; this implies that, on the 
bottom triangle, one has 




D\h + Dlh 


H^A^B^) 

TPpoipy ’ 


2 


D\h + D^h 


H^{A^A^) 

pi (AMI) ’ 


where by D\h and D\h we denote the constant value of Dih and D 2 h on the bottom 
triangle. This readily implies 


^ U^[A^B^) ^ pi(AMi) “ f 


( 2 . 6 ) 


On the other hand, 

V}{A^B^)= \Dg\<m} , 


H\A^A^) = \Dg\dH\ 
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which inserted in (2.6|), and using ( 2 . 2 ) from Step 1, gives 


\D’^h\ < - 


r JB{Vi,f)ndSo 


iDgldn^^K iDgldn^. 


'dSo 


Hence, we deduce that 

\Dh\ = 






D^h\ < 


Kf 


^t 2 


'dSo 


\Dg\d'H^<K \Dg\dV}, 


'dSo 


thus the validity of (2.5) follows. 

Let us now consider the degenerate case, where in the bottom triangle the function h 
is made by two degenerate affine pieces, one on the left triangle A^PA^ and the other on 
the right triangle A^PB^] we call and D^h the constant values of Dh respectively on 
the left and on the right triangle. Since the image of the segment A^B^ through the map h 
is the path 7 ^ (that is, the union of the two segments and A^B^), parametrized at 

constant speed, we get that \D[h\ = \D[h\ (while in general D[h 7 ^ D\h)] more precisely. 


\D{h\ = \Dlh\ = 


V}{A^A^) +V}{A'^B^ 


(2.7) 


Moreover, the affine map in the right triangle transforms the segment A^B^ in the segment 
A^B^^ while the affine map in the left triangle moves A^A^ onto APA^] this implies 


V2 


Dill + Dili 


hHa^b"^ 


n^[A^B^) ’ 


2 


D[h + Dlh 




P^{A^A^) ’ 


which together with (2.7) gives 
3 r 


\DP\ < 


r JB(Vi,f)ndSQ 


\Dg\dH\ 


\D^h\ < 


r JB{yi,f)r\dSQ 


\Dg\dP^ 


Arguing exactly as before, again thanks to ( 2 . 2 ) of Step 1, we obtain again the validity 
of (2.5), possibly with a slightly larger, but still purely geometric, constant K. 

Step 9. Estimate for !vADh\- 

In this step, we want again to find an estimate for the integral of \Dh\, but this time on 
the generic horizontal strip "Dj, 1 < i < k — 2. Our goal is to obtain the estimate 


[ \Dh\<K\Vi\[ \Dg\dH^ + K[ \Dg\dn\ 

Jv, JSo 7a»A'+iuB»S*+i 


(2.8) 


where by \Vi\ we denote the area of the horizontal strip Consider the horizontal 
segment by symmetry, it is not restrictive to assume that this segment is below 

the x-axis, precisely at a distance 0 < r < 1 from the “south pole” Vi = ( 0 , — 1 ); in other 
words, we have that = (—r, r — 1) and = (r, r — 1). Moreover, let us call a the 

distance between the segment and the segment A^B^, and 


i:=max{n^{A^A^+^),n\B^B^^^)] < [ \Dg\dn\ (2.9) 


Remember now that, in Step 7, we defined h as the function which is affine on each of the 
triangles PjPj+iQj, and Pj^iQjQj+i, sending each of the points point P^ (resp., Qm) in 
Sq onto Pm (resp., Qm) in «So. Let us then concentrate ourselves on the generic triangle 
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T = PjPj+iQj (for the triangles of the form Pj+iQjQj+i the very same argument will 
work); since h is affine on T, let us for simplicity denote by D'^h the constant value of 
Dh on T. 

First of all let us recall that, on the segment the function h has been defined 

as the piecewise linear function whose image is parametrized at constant speed; this 

ensures that 


Dlh 




( 2 . 10 ) 


We observe now that, by definition, is the shortest path in the closure of <So con¬ 
necting with in particular, is shorter than the image, through g, of the 

curve connecting with on dSo passing through the south pole. This implies in 
particular that 




lB(Vi,V2r)naSo 


iDgldn^ 


which inserted in (2.10) and recalling (2.2) gives 


\Dlh\ < 


-f 

2r JB(Vi,V2r)ndSo 


\Dg\dn^<K \Dg\dn^ 


'So 


( 2 . 11 ) 


Let us now use the fact that the affine map h on T sends the segment PjQj onto the 
segment PjQj- Calling, as in Figure]^ d and d' the distances between and Pj, and 
between A* and Qj, we derive that 


[d! -d + a)Dlh + a Dlh = H\PjQj) < 


( 2 . 12 ) 


where in the last equality we have used the estimate (2.4) from Step 6, which is valid 
since PjQj is a vertical segment in the sense of Step 6. 



^i+i 





Figure 4. Position of points and lengths in Step 9. 


Let us now use once again the fact that is the shortest path between and 

pi+l 

on the closure of «So: in particular, y®+^ is shorter than the path obtained as the 
union of A*’''^A*, the part of y* between A* and Qj, the segment QjPj-, and the part of 
between Pj and namely. 


pi(y*+^) <-Hi(A*+'A*) + 
d'V}{-i^) 


d'V}{-i^ 




+ V}{Q^P,)+V}{^ 


d 




-IV 


^ / 1 


d 
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which implies 


d 


W(Y») ,, 


The completely symmetric argument, using that 7 * is the shortest path between A* and 
B\ thus shorter than the union of the part of between A*"*"^ and Pj, the 

segment PjQj, and the part of 7 * between and B\ gives the opposite inequality, hence 
we get 


d 


-d' 




< 2i, 


which further implies 


|d-d'| 


n\A^B 


<2i + d 


) 




<2i + 
<2i + 


- n\'y 


•^ 1(^+1 




+ 2a 


Using now (2.10), we can rewrite the above estimate as 

\d-d'\\Dlh\ < Ai + 2a\Dlh \, 


< 4£ + 2cr 


W(7 


W‘(yf+‘B‘+>) ' 


which recalling also ( 2 . 12 ) finally gives 

o’lDlh] < 5£ + 3a\Dlh\ 


We can then easily evaluate the integral of |Zlh| on T, also by (2.11), as 


1^\Dh\ = \D^h\ < \Dlh\ + \Dlh\ < in 


AK 


/ n^|d-H^ + 5- 

'So O' 


Adding now the above estimate over all the triangles T forming and recalling (2.9) 
and the fact that \Vi\ < a, we directly obtain ( 2 . 8 ). 


Step 10. Definition of the modified function h and conclusion of the proof. 


We start observing that, adding the estimates (2.5) for the top and for the bottom triangle 


tog ethe r with the estimates ( 2 . 8 ) for all the horizontal strips, we directly obtain the va lidit y 
of (2.1) for the function h. However, the proof is still not over, because h satisfies ^2.1), 


coincides with g on dSo, and it is hnitely piecewise affine, but it is not a homeomorphism 
(unless all the paths 7 ® lie in the interior of «So). However, we can easily obtain this with 
a simple modification of h: more precisely, let us slightly modify all the paths 7 ®, so that 
they remain piecewise linear but they live in the interior of <So and they do not intersect 
each other. The idea, depicted in Figure is obvious. Notice that, since there are only 
hnitely many paths 7 ®, and each of them has only hnitely many vertices, it is clear that we 
can “separate” as desired all the paths, and we can also move each of them of a distance 
which is arbitrarily smaller than all the other distances between extreme points. Then, 
we dehne the function h exactly in the same way as we dehned h, except that we use 
not the original paths 7 ® but the modihed ones; therefore, the function h is now not only 
hnitely piecewise affine and coinciding with g on OSq, but it is also a homeomorphism. 
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Figure 5. Modification of the paths 7* in Step 10. 


Moreover, the estimate (2.1) is still valid, with a geometric constant K which is as close 
as we wish to the one found above. The proof of the theorem is then now concluded. □ 

Remark 2.2. A trivial rotation and dilation argument proves the following generalization 
of Theorem 2J_. If S is a sguare of side 2r and g : dS —)■ is a piecewise linear and 

one-to-one function, there exists a piecewise affine extension h : S ^ of g such that 

( \Dh{x)\dx<Kr ( \Dg{t)\dH\t). (2.13) 

JS JdS 


3. Extension in the degenerate gase \Df{c) \ = 0 but J/(c) 7^ 0 


As already explained in the description of Section fM, a crucial difficulty in our proof 
will be the case when a square S is “good” (this means that Df is almost constantly equal 
to some matrix M within S), but detM = 0, while being M 7^ 0. It will be important 
to handle this case with care, because the map / on 5 is then very close to an affine 
map, but this affine map is degenerate. The goal of this section is to prove a single result, 
which will solve this difficulty. Recall that, whenever a map g is dehned on dS, for any 
t G dS we denote by r(f) the tangent vector to dS at t, by Dg(t) the derivative of g in 
the direction r(t) (whenever it exists), and by Ids \Dg \ dTL^ the length of the curve g{dS). 


Theorem 3.1. Let S be a sguare of unit side and g : dS a piecewise linear and 

one-to-one function such that 



■ 


dV}{t) < 6 


(3.1) 


for 6 < (^max? where 5max <C 1 zs a geometric guantity. Then there is a finitely piecewise 
affine homeomorphism h : 5 —)■ such that h = g on dS and 



dx<K6. 


(3.2) 


Proof. We divide this proof in several steps, to make it as clear as possible. 

Step 1. Definition of good and had intervals. 

First of all we notice that, since g is a. one-to-one piecewise linear function, then its image 
g{dS) is the boundary of a nondegenerate polygon, that we call S. Moreover, thanks 
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to (3.1), we know that this polygon is very close to a horizontal segment in Up to 


a translation, we can assume that the hrst coordinates of the points in S are between 0 
and L. 

Fix now any 0 < a < L: it is reasonable to expect that there are exactly two points 
in g{dS) having hrst coordinate a, and that the two counterimages in dS are more or 
less one above the other (this means, with the same hrst coordinate). In the situation of 
Figurethis happens with a, but not with a', since the points in dS with hrst coordinate 
a' are four. We dehne then “good” any a G (0, L) with the property that exactly two 
points in dS have hrst coordinate a, and with the additional requirement that a > 26 
and a < L — 25. For any such a, we call and the two above-mentioned points. 



9{dS) 



Qa\ 

— - 


o' 

a 

a' 

L 


Figure 6 . A good a and a bad a' in Step 1. 


being above and we call = g ^{Pa) and Qa = g ^{Qa)- can immediately 
show that a big percentage of the points are good, more precisely 


W 


jcr G (0, L) : a is not goodj^ < 56. 


(3.3) 


Indeed, take any segment RS in dS on which g is linear, and call as usual R = g{R) and 
S = g{S): by dehnition, we have that 


IRS 


J Q ) ■ 


dR}{t) > 
> 


/ Dg{t)dH\t)- 

JRS JrS 

(5i-Pi)-(Fi-Pi)| 


nit) 

0 


dn\t) 


(3.4) 


As a consequence, if the segment RS is going backward (that is. Si — Ri and Si — Ri 
have opposite sign), then its horizontal spread is bounded by the above integral on the 


interval RS. Recalling (3.1), this means that all the backward segments have a projection 


on (0, L) with total lenght less than 6. Since of course any a G (25, L — 26) which is not 


good must belong to this projection, the validity of (3.3) follows 


By adding the inequality (3.4) for all the segments of dS, we hnd also that L equals 


the horizontal width of S up to an error 5/2, thus in particular 

5 ^5 

1- - < L < V2+ - . 

2 ~ ~ 2 

Moreover, take any good a and consider all the segments on dS connecting P„ and Q„: 


again adding (3.4) on all these segments, and recalling that Po- and have the same 
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first projection, we derive that 


{Pa)l-{QM < 


2 ’ 


(3.5) 


that is, the points Po- and Qo- are always exactly one above the other np to an error 5/2: 
the factor 1/2 comes by the possibility of choosing either of the two paths in S from P^ 
to (j. 

Finally, assnme that P^ and lie on a same side of S for some good a. Adding once 


again (3.4) among all the segments where g is linear connecting and Qo-, we derive that. 


up to an error 5, the sum of all the horizontal spreads of these segments coincides with 
the corresponding sum of the horizontal spreads in dS] however, the hrst sum is smaller 


than 5/2 by (3.5), while the second is at least the minimum between 2cr and 2(L — a) 


which is impossible by the definition of good a. In other words, we have proved that Pa 
and Qa never he on a same side of 5 if cr is good. 

Observe now that, since g is piecewise linear, then by definition (0, L) is a finite union 
of intervals, alternately done entirely by bad a and entirely by good ones. However, the 
endpoints of all these intervals are always bad. Therefore, we slightly shrink the intervals 
made by good points and we call good intervals these shrinked intervals. Thanks to ( 3.3[ ), 
we can do this in such a way that the union of the good intervals covers the whole (0, L) 
up to a length of 65: notice that all the points of any good interval are good points, 
also the endpoints, while the bad intervals may also contain good points. Finally, it is 
convenient to make the following further slight modification: up to replace a good interval 


with a finite union of good intervals, we can also assume that whenever 


cr,a 


is a good 


interval, the map g is linear in the segments PaPa' and QaQa'- 

Step 2. The extension on the segments PaQa, definition of good and bad guadrilaterals. 
In this step, we extend g -which is defined on dS- to a union of segments in S. More 
precisely, recall that (0, L) has been divided in intervals, which can be either bad or good. 
Moreover, the extremes of these intervals are all good points except 0 and L. Take then 
any other of these extremes, say a, and consider the points Pa and Qa in dS. We define g 
on the segment PaQa as the linear function such that g{Pa) = Pa and g{Qa) = Qa'^ notice 
that all the different open segments PaQa are contained in the interior of S by Step 1, and 
they do not intersect with each other by construction. We have then many segments PaQa 


inside S, almost vertical by (3.5), on each of which a linear function g is dehned. Observe 


that, as a consequence, S has been divided in several quadrilaterals (actually, the hrst 
and the last one are generally triangles), and g is dehned in the whole corresponding 1- 
dimensional grid; also S has then been subdivided by the images of in a union of several 
polygons. A positive consequence of this fact is that we can now dehne the extension h of 
g in an independent way from each quadrilateral in S to the corresponding polygon in S 
-respecting of course the boundary data, and being a piecewise affine homeomorphism: 
then, the resulting function h will automatically be a piecewise affine homeomorphism. 

Let us conclude this short step with another piece of notation: any quadrilateral in 
S will be called a good guadrilateral if it corresponds to a good interval in (0 ,L), and a 
bad guadrilateral otherwise. In the remaining of the proof, we will hrst give an estimate 
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for the good quadrilaterals; then, we will give one for the hrst and the last quadrilateral, 
that is, the one starting at 0 and the one ending at L: notice that these quadrilaterals 
are always bad by dehnition, and actually they are usually triangles. Finally, we will give 
an estimate for the “internal” bad quadrilaterals, which we obtain by considering two 
subcases. 

Step 3. The extension in good quadrilaterals. 

Let us hrst start by considering a good quadrilateral, corresponding to the good interval 
(cr, ex') in (0, L); for brevity, we will write P, P', P, P' in place of P^, Pa', Pa, Pa'- Recall 
that the map g is linear between P and P', as well as between Q and Q', thanks to the 
construction in Step 1. As a consequence, the image under g of the boundary of the 
quadrilateral PP'Q'Q is the boundary of the quadrilateral PP'Q'Q, and we have to 
dehne the extension h by sending the interior of PP'Q'Q onto the interior of PP'Q'Q. 
Let h simply be the piecewise affine map sending PP'Q onto PP'Q and QP'Q' onto 
QP'Q'. 



Figure 7. The approximation in a good quadrilateral. Step 3. 
We need to estimate 

r /I n\ 

dx , 


IPP'Q 


Dh- 


0 0 


(3.6) 


the estimate in the triangle QP'Q' being then of course identical. Let us dehne for 
shortness 

i = P 2 -Q 2 , b = P[-P,, ^ = P' - P 2 , 9 = arctan(e/6), 

6i = Pi-Qu a = P 2 -Q 2 , v = P',-Pi, /3 = P'^-P 2 , 

we refer to Figure for help with the notation. By dehnition, the constant value of Dh 
in PP'Q satishes 

bDihi + ^D2hi = 7] , 
diDihi + iD2hi = 0 , 

Let us start by dehning 


bDih2 + ^D2h2 — fd , 

6iDih2 + £^ 2^-2 = Oi . 


(3.7) 


£ = 


IPP' 


Dg{t) - ( J 0 ) ■ 




(3.8) 
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SO that adding the values of e on the different segments we will get less than 5 by (3.1). We 


claim now the validity of the following estimates, all obtained again arguing as in (3.4): 

A 

\v 




< e, a < S, 


5 maxjtan^, 1} . (3.9) 


The hrst two estimates can be found just integrating on the segment PP', so they are 
valid with the small constant S] instead, to get the third estimate we have to integrate on 
all the segments connecting P and Q on dS, so we can only estimate with 6] the fourth 


estimate is given by (3.5). Finally, the evaluation of i follows by a simple geometric 


argument, just recalling that a > 2S, (3.5) and that we have dehned 9 as the direction of 
the side containing PP'. 


Let us now start by evaluating Dihi\ inserting the third equation of (3.7) into the 
hrst one, we get 


DM b 


= h: 


from which we readily obtain, by using the estimates (3.9) and recalling that ^/b = tan6', 

(3.10) 


\Dihi 


I ^ 

<21 - + — 
- ' b bi 


Substituting the value of Dihi again in the third equation of (3.7), we get then 


£> 2^1 


i 


Dih^ 


6 e ^6 

< -^-h — 

- 2i b bi 


(3.11) 


We control now the derivatives of h 2 '- inserting the second equation of (3.7) into the 
fourth, we get 


712^2 


.a M 


SO that, again using ( 3.9[ ) and again recalling that ^/b = tan^, we deduce 
7^2 ^2 


5 5e 

< 2 - H-, 

- 7 67 ’ 


Dih-} 


e ^6 
< 2 - + 2 — . 
- b be 


(3.12) 


Estimating the integral in (3.6) is then straightforward. Notice that 5 is a hxed constant. 


not depending on the subdivision in intervals: as a consequence, we can assume without 
loss of generality that ^ < S < i, otherwise it is enough to subdivide a good interval in a 
hnite union of good intervals; the area of the triangle PP'Q is then less than 67, and so 


from (3.10), (3.11) and (3.12) we obtain 


tPP'Q 


Dh 


(oo) *s(5^ + 5| + 3^ + |)-«<8<J + y5? + 36). 


where we have also used that 6 < 1/2 and 7 < 3/2 (the latter follows by straightforward 
geometrical arguments). Of course, the fully analogous estimate holds for the integral in 


the triangle QP'Q', up to replace the segment PP' by QQ' in the dehnition (3.8) of s. 


To conclude, we need to evaluate the total integral in the union of the good quadri¬ 
laterals; this is simply achieved by summing the above estimates over all the different 
quadrilaterals. Notice that the constant 6 is hxed and does not depend on the quadrilat¬ 


eral, while the constants e, 6 and ^ are specihc of each quadrilateral. By dehnition (3.8) 
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of e, it appears clear that the sum of all the different e’s is less than 6, while by dehnition 
of the lengths in the square it is clear that the sum of the different as well as of the 
different b, is bounded by 4. As a consequence, we deduce that 



dx < K6 , 


(3.13) 


where G denotes the union of all the good quadrilaterals in S, while K is as usual a purely 
geometric constant. 

Step 4- The extension in the first and last bad quadrilaterals. 

In this and in the next step we are going to consider the bad quadrilaterals. Notice that, 
since almost the whole square is done by good quadrilaterals thanks to (3.3), we can even 
be satisfied by a rough estimate here, while we needed a precise one in the preceding step: 
what is important, is that we can dehne a piecewise affine homeormophism h on each of 
the bad quadrilaterals. 

Here we begin with the “first” and the “last” quadrilateral, that is, with the quadri¬ 
laterals which correspond to the two intervals having 0 or L as one endpoint. Notice that 
these “quadrilaterals” are actually triangles, unless some side of the square is very close 
to being vertical. More precisely, let us consider just the first bad quadrilateral C, by 
symmetry: as Figure [^depicts, it can be either a triangle VPQ, being V the left vertex of 
the square, or a quadrilateral VV'PQ, if V and V are the two left vertices of the square, 
being the side W almost vertical, and then the sides VP and VQ almost horizontal. 
Notice that all the sides of C belong to dS except PQ. We need to send C on the polygon 




Figure 8 . The approximation in the hrst bad quadrilateral. Step 4. 


C inside S made by the points which have first coordinate less than a = Pi = Qi, shaded 
in the right of the figure. Keep in mind that by construction (recall Step 1) the coordinate 
a is good, and the bad intervals cover only a portion less than 25 of {26, L — 26): this 
means that 26 < a < 46. As a consequence, again by using several times (3.4) and (3.5), 
we know that 




2 <^1-^1 <55, 


^<Qi-Vi<56, 


\Qi-Pi\<^, 


and the estimates on V are valid also for V in the case when the bad quadrilateral C 
is actually a quadrilateral. We would like to infer that C is the biLipschitz image of a 
square with side 6, with uniformly bounded biLipschitz constant; however, this is true 
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onl y if i is comparable to 6, while we only know by Step 3 that £ > 6 -this was established 
in (3.9). Let us then consider the affine map ^{xi^X 2 ) = (xi,X 2 ( 5 /^), and let C = ‘h(C); 
dehne also g = g o on dC, which is admissible since g is dehned on the whole dC. 
By construction, the set C is the biLipschitz image of a square of side 6, with biLipschitz 
constant less than a geometrical constant K. We can then apply Theorem |2.1| to the map 


g, and we hnd an extension h of g inside C such that ( 2 . 1 ) holds, that is. 


hDh{y)\dy<K6 l\D~g{t)\dn\t) 

Ic JdC 


(notice that the multiplication by 6 comes from the argument of Remark 2.2). Observe 
now that the integral in the right side of the above inequality is simply the perimeter of 


C, which is less than K5 by (3.1) and again by (3.4). Thus, we infer 


\Dh{y)\dy < KS^ 


(3.14) 


Finally, we dehne h = h o 4) on C: this is a piecewise affine homeomorphism from C to 
C, and by dehnition it extends the map h already dehned on dC. We have then to show 
that Dh is not too big on C, and to do so it is enough to observe that 

\Dh{<!>-\y))\ < \Dh{y)\, 


which by (3.14) hnally implies 


\Dh{x)\dx < J^\Dh{^-\y))\^dy <‘^ J^\Dh{y)\dy < K6 . (3.15) 

We have then found the estimate we were looking for related to the hrst bad quadrilateral, 
and by symmetry the same holds also in the last bad quadrilateral. 

Step 5. The extension in the internal bad quadrilaterals. 

To conclude our analysis, we need to concentrate in the internal bad quadrilaterals. Let 
C be a bad quadrilateral, and let us call its vertices, as usual, P, Q, P' and Q'; the image 
of dC under g is then the boundary of a polygon C. Notice that C needs not to be 
a quadrilateral, since it has two vertical sides, PQ and P'Q\ but PP' and QQ' are 
piecewise linear paths, not necessarily segments. Keeping a notation similar to that of 
Step 3, we set 

£ = P2 - Q2 ■, b = P[- Pi, a = max{P2, P'2} - Lnin{Q2) Q2} ) 

^ = P'-P 2 , 0 = arctan(^/6), g = H\PP') + , 

see Figures and 10 By a simple symmetry argument, we can assume without loss of 

generality that 

6 > 0, and either PP' and QQ' are parallel, or 0 > 7r/4. (3.16) 

Observe that this is possible because, if PP' and QQ' are not parallel, then they belong 
to two consecutive sides of the square, hence if 0 < 7r/4 we just have to exchange P with 
Q. Notice that, by dehnition. 


V = 


IPP'UQQ' 


Dg(t) dn\t). 


(3.17) 



















DIFFEOMORPHIC APPROXIMATION OF PLANAR SOBOLEV HOMEOMORPHISMS 25 


We need now to further subdivide our analysis in two subcases, depending whether a is 
bigger or smaller than lOr]. Notice that a is bounded by 6, while r] could be even much 
smaller than 6, since the sum of all the different r/’s corresponding to bad intervals is 
smaller than 35: indeed, the total length of the internal bad intervals is less than 25, so 

1 O' 


we do not even need to subtract the matrix (^g g J as in (|3.8|). As a consequence, either 

of the two cases may actually hold. 

Step 5a. The case a < lOr/. 

Let us start with the case when a < lOr;. We let H be the point in the segment P’Q' 
satisfying P 2 = PI 2 (such a point exists by (3.16)), and H = g{H), which is well dehned 
since g has been dehned on the good segment P'Q'. We subdivide the quadrilateral C 



Figure 9. The approximation in an internal bad quadrilateral: case 1, Step 5a. 


into the union of the triangle PP'H and the quadrilateral PHQ'Q, and we aim to dehne 
the function h separately one these two pieces. First of all, similarly as in the proof of 
Theorem 2T, we consider the shortest path between P and H in C, which is a piecewise 
affine path, possibly intersecting dC in other points than P and H, and we call 7 a 
slight modihcation of this path, which is still piecewise affine, but which is entirely in the 
interior of C except for the two extremes P and H. By minimality, we can of course take 
the modihed 7 satisfying 


W‘( 7 ) < n\PP') + H'-(P’H 


(3.18) 


We extend then g to the segment PH as the piecewise affine function sending the segment 
PH onto the path 7 at constant speed. 

Let us now point our attention on the triangle PP'H\ the segment PH is horizontal 
by dehnition, while the segment P'H is “quite vertical”; more precisely, it is contained in 
the segment P'Q' and by dehnition we have 

\P[-Q[\<1, P'-Q',>i>6. 

The triangle would then be a biLipschitz image of a square with side b, with uniformly 
bounded constant, if ^ were not too much bigger than b, or, in other words, if 6 is not 
too big. Since we cannot be sure that this is the case, exactly as in Step 4 we dehne <h 
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the affine map which does not modify the horizontal segments, and which shrinks of a 
ratio ^/h the segments parallel to P’H. Then, ^{PP'H) is a triangle which is nniformly 
biLipschitz with a sqnare of side b, so exactly as in Step 4 we apply Theorem im to the 
map g = go hnding an extension h on ^{PP'H), and we finally obtain the extension 
to g in PP'H as h = h o $. Estimating the derivatives of h, h, g and g exactly as in 
Step 4, we get then the estimate 

f \Dh{x)\dx < f \Dh{y)\dy < f \Dg{t)\dt 

JPP'H b J^{PP'H) Jd{<S>{PP'H)) 


= K^V}[d{g{PP'H))) = Ki[v}{PP')+V}{-i)+V}{P'H 


< K['H\PP')+'H\P'H)) < K{ri + a) < Ky, 


where we have also nsed (3.18) and the assnmption a < lOy. 

Let ns now consider the qnadrilateral PHQ'Q. Since we have already seen that 
PQ and P[Q' are “qnite vertical”, while PH is exactly horizontal and QQ' is “qnite 
horizontal” since it makes with the horizontal direction an angle eqnal to n/2 — 6 < ir/A, 
this qnadrilateral is nniformly biLipschitz with a rectangle. Up to shrinking vertically of 
a ratio b/i as before, it becomes nniformly biLipschitz with a sqnare of side b, hence by 


argning as before by shrinking, applying Theorem |2.1| and then stretching back, we define 
an extension h of g inside the qnadrilateral PHQ'Q which satisfies 


/ \Dh{x)\dx < KiiH^PQ) +H\HQ') +H\y) +H^QQ')) < Ky, 

JPHQ'Q ^ '' ' J 

which pnt together with the estimate above for the triangle PP'H gives 


\Dh{x) \ dx < Ky . 


( 3 . 19 ) 


Step 5h. The case a > lf)y. 

The last case that we have to consider is that of an internal bad qnadrilateral correspond¬ 
ing to a > lOr^: in this case, if we argned as in Step 5a, then we wonld find the same 


estimates as in (3.19), bnt with K5 in place of Ky] and in tnrn, this wonld not be ac¬ 


ceptable, becanse adding all the different r/’s for the bad qnadrilaterals we get something 
smaller than 5, while adding a term 5 in each of the bad qnadrilaterals we conld get any 
big constant in the end, since bad intervals conld be many more than 1/5. 

As a conseqnence, in this snbstep we present a different dehnition of the extension h for 
the case a > IDy. More precisely, take fonr points iTi, iT 2 , as in Fignre [T0| in 

the segments PQ and P'Q', at a distance y from the fonr vertices; let also Hi = g~^{Hi) 
for i = 1, 2, 3, 4. By definition of y, the open segments H 1 H 2 and H 3 H 4 are entirely 
contained in the interior of C, and by constrnction the same happens for the segments 
H 1 H 2 and H 3 H 4 in C. We regard then both C and C as the nnion of three pieces: 
the internal qnadrilaterals HiH 2 H 4 ^H^ and H 1 H 2 H 4 H 3 , and the “top” and “bottom” 
remaining parts, shaded in Fignre We aim to define the piecewise affine fnnction h so 
to send each part of C onto the corresponding one in C. 
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Figure 10. The approximation in an internal bad quadrilateral: case 2, Step 5b. 


For the “top” and “bottom” part, we can argue more or less exactly as in the last 
steps: each of the quadrilateral PP'H 2 Hi and H^H^Q'Q can be transformed into a square, 
then one applies Theorem 2.1| and then goes back to the quadrilateral; since the perimeter 
of each of the shaded regions in C is now at most Arj, the same estimates as in Step 5a 


can be repeated, so that similarly to (3.19) we get now 


/ \Dh{x)\dx < Kr]. (3.20) 

JpP'H2HiUH3HiQ'Q 

To conclude, we have to dehne the extension h so to send the internal quadrilateral of C 
onto the internal quadrilateral of C, and we will do that again by sending in an affine way 
the triangle H 1 H 2 H 3 (resp. H 3 H 2 H 4 ) onto the triangle H 1 H 2 H 3 (resp. H 3 H 2 H 4 ). We 
need thus only to check the value of \Dh\ on the triangle H 1 H 2 H 3 , being the estimate for 


the triangle H 3 H 2 H 4 completely similar. To get the estimate, as in Figure 10 we set 


H4-H3 = {Sui), 


H2-H4 = (b,0, 


9 = arctan C 
b 


notice that (5i, £, b, ^ and 6 are very similar to (5i, £, b, ^ and 6 , since the assumption 
a > lOrj implies that the points Hi are very close to the vertices of C. Hence, the constant 
matrix Dh in H 1 H 2 H 3 satishes 


(a — 2 ? 7 , 0) < S , 


\Dh{ 6 iJ) 

|Z)h(6, 6tan0)| = |Zlh(&,^) 
As a consequence, we get hrst that 


= \H2 


Hi\ <r]. 


(3.21) 


c 

Dh{Si, tan 0)1 < rj < , 

' b b 


and then that 


\Dh{0,£ 


(5itan0)| <5(1 + ^ 


(3.22) 


Recall now that the estimates (3.9) ensure 


^ = — > 2 tan 6 > 2 tan 9 ; 
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notice carefully that the estimates (3.9) were obtained in a good quadrilateral, so they 


are not valid now, but since the segment PQ corresponds to a good a, and in particular 
it is in the boundary of the good quadrilateral immediately preceding C, the estimates 


about i and 5i are still valid and we can use them now. As a consequence, by (3.22) we 
get 




(3.23) 


and substituting this in (3.21) we have also 


Dh{b,0) <ri + btan 6 \D 2 h\ < rj + 2btan6 ^ j < rj + 36(l + ^ ) < 4?7 + 36, 


from which we derive 


iDihl < 5 ? + 3 

I I 6 


(3.24) 


Since the area of the triangle H 1 H 2 H 3 is bounded by b{i + btanO), by (3.23) and ( |3.24[ ), 
and recalling again 6 tan 6 < i, we get 




\Dh(x) \ dx < (bi + 6^ tan 0 ) (^3 ^ + 3^ + 5^ + 

^ \ i bi b J 

< 35b + 36^ + 357] + 3r]b + 3r]i + 5r]^ + 3bi + 36^ < K{b + t]) , 


so that repeating the same estimate in the triangle H 3 H 2 H 4 , and adding (3.20), we obtain 
that in a bad quadrilateral C where a > lOr] the estimate 


\Dh{x) \ dx < Kb + Krj 


(3.25) 


holds. 

Step 6. Conclusion. 

We can now put together all the estimates of the last steps to conclude the proof. Let 
us start by considering the bad quadrilaterals, whose union is 5 \ G, since in Step 3 we 


have defined G as the union of the good quadrilaterals. Thanks to (3.19) and (3.25), we 


know that the integral of \Dh\ in any internal bad quadrilateral can always be estimated 
by 6 + r;. If we add the different 6’s corresponding to the bad quadrilaterals, up to an 
error 5 we find the sum of the lengths of the internal bad intervals, which is at most 


25 by construction. On the other hand, adding the different p’s and recalling (3.17), we 


get something smaller than (twice) the sum of the lengths of the bad intervals. As a 
consequence, putting together the estimates for all the internal bad quadrilaterals, and 


also adding the estimate (3.15) for the first and the last bad quadrilateral, we obtain 


/S\G 


inh(x)ldx < K5. 


Since the total area of the bad quadrilaterals can be estimated by (twice) the total length 
of their horizontal projections, which in turn corresponds with the total length of the bad 
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intervals up to an error 6, and so it is less than 76, we can now insert (|3.13|) to get 
1 0 " 


Dh{x) - 


0 0 


dx = 


JG 

<K6 + 


Dh{x) — 


1 0 
0 0 


dx + 


ls\G 


Dh{x) - 


1 0 
0 0 


dx 


[ \Dh{x)\dx+ S\G <K6, 
Js\G 


which is (3.2), and then the proof is concluded. 


□ 


Remark 3.2. A trivial rotation and dilation argument proves the following generalization 
of Theorem 3J_: whenever S is a sguare of side 2r, g : dS is a piecewise linear 

and one-to-one function, and M is a matrix with det M = 0, there is a finitely piecewise 
affine extension h : S ^ of g such that 

[ Dh{x) — M dx < Kr [ Dg{f) — M ■ rif) dVfif), 

Js Jas 

as soon as 

f \Dg{t) - M ■ T{t)\dV}{t) < r(5MAx||M|| . 


4. Proof of Theorem 11.11 


This last section is devoted to give the proof of Theorem 1^ This proof is still quite 
involved, but the overall idea is simple after the preceding sections. As already explained 
in the introduction, the idea is to divide the whole fl in squares, and then treat them in 
three different ways: roughly speaking, the “good” squares, where the function is very 
close to an affine map (and this group will be further divided in two subgroups), and the 
“bad” ones, where this is not true. Moreover, we will have to slightly change the value 
of / on the boundaries of all these squares, in order to become piecewise linear. Then, in 


the bad squares we will simply use Theorem 2.1 to get an extension, and the constant K 


in (2.1) will not be a problem because the bad squares will cover only a small portion of 


fl. Instead, we have to perform a very precise approximation of / in the good squares; 
to do so, we will treat differently the squares where the affine map close to / has zero 
determinant, and those where the determinant is strictly positive. For the hrst ones, we 


will use Theorem |3.1[ while for the second ones it will be enough to interpolate the values 
of / on the boundary, as we show in Section 4.1[ 

Before starting with the proof, let us give a couple of dehnitions. 


Definition 4.1. We say that S{c,r) is a Lebesgue square with matrix M G and 
constant 5 > 0 if S{c,3r) C fl and 


4 - \Df{z)-M\ 

Js{c,3r) 


dz < 6 . 


Definition 4.2. Let S{x,r) O Q be a sguare, and denote by Ti and T 2 the two triangles 
on which S is subdivided by the diagonal connecting (xi — r,X 2 r) and (xi -\- r,X 2 — r). 
We call g:>s(x,r) the piecewise affne function which is affne on the two triangles Ti and 
T 2 , and which coincides with f on the four vertices of S{x,r). 
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4.1. The Lebesgue squares. In this first subsection, we consider the situation in the 
best possible case, namely, a Lebesgue square. It is rather easy to show the following 
uniform estimate. 


Lemma 4.3. For every e > 0 and every matrix M, there exists 6 = 5{M,e) <C e such 
that the following holds: if S{c, r) is a Lebesgue sguare with matrix M and constant S < 6, 

re 


11/ - V?||L“(5(c,r)) < re , \\Df - D(f\\Li(Sic,r)) < r^e , ||/ - 'lp\\ L^(S{c,2 r)) < 


10 


-)■ 


(4.1) 


IS an 


where ip = ips(c,r) is the piecewise affine map of Definition f.2 and : 
affine function satisfying Dfi = M. If in addition det M > 0, then ip is injective and 

- e)r)) Cip{S{c,r)) C f[S{c,{l + e)r)). (4.2) 

Proof. We assume for simplicity of notation that the point c is the origin of and we 
write iS(r) instead of 5(0, r). Let i? be a big constant, depending only on M and e, to be 
specihed later, and let us dehne the two sets A and B as 

A = |x e (-3r, 3r) : -P" \Df{x, t) -M\dt> , 


B = 


J-3r 

r3r 


{y e (-3r, 3r) : \Dfif y) -M\dt> . 


By dehnition of Lebesgue square, we immediately get that 

6r 


1^1 < 


R 


B < 


6r 

1 r 


(4.3) 


Let us now arbitrarily fix a point z = {x,y) E S(r) with x ^ A, y ^ B, and let us dehne 
: S{c, 3r) —)■ as = f{z) + M{w — z): it is clear that ip is an affine map with 

Dip = M, and calling g = f — ip we have by dehnition g{z) = 0. We claim that 

\g{w)\ < 12rR5 \/w = {x,y) E 5(3r) \{Ax B). (4.4) 

Indeed, assume that x ^ A {if y ^ B the obvious modihcation of the argument works). 
Recalling that g{x, y) = 0 and that x ^ A and y ^ B, we can evaluate 

\g{w)\ = \g{x,y) - g{x,y)\ < \g{x,y) - g{x,y)\ + \g{x,y) - g{x,y)\ 

< r\Df{x, t)-M\dt+ r \Df{f y) -M\dt< l2rR5 , 

Jv J X 


and (4.4) is proved. 

Let now w = {x.,y) E S(2r) be a generic point. By ( 4.3[ ), we can hnd xi < x < X 2 
and yi < y < y2 such that for i = 1, 2 we have {xi,yi) E 5(3r), Xi ^ A, yi ^ B and 
X 2 — xi < 6r/R, 1/2 — 1/1 < 6 r/R. Hence, w is inside the small rectangle TZ having sides 
with coordinate Xi and yi, and by (4.4) we know that 

\g{P)\<l2rR5 y P E dlZ. (4.5) 

By dehnition ip{dP) is a small parallelogram around ip{w), with 

ip{dn) C b(^iP{w), ||M|| 
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having defined ||M|| = max |M(t)|/|t|. This estimate, together with ( 4.5[ ), ensures that 
the whole curve f{dTZ) is done by points with distance less than 12rR5 + 6r\/2\\M\\/R 
from ipiw). Since / is a homeomorphism, the point f{w) is inside this curve, hence we 
finally deduce 


|/(to) - = | 9 (m)| < 12rRS + pf|| < r A _ 


( 4 . 6 ) 


where the last inequality holds as soon as R has been chosen big enough, depending on 
M and on e, and then 5 has been chosen small enough, depending on R and e, and thus 


ultimately only on M and Hence, we have obtained the third estimate in (4.1). 

Let us now consider the function 93 = (fsir), and let us call 

Hi = (ci - r, C 2 + r), V 2 = (ci + r, C 2 + r), 1/3 = (ci + r, C 2 - r) , V 4 = (ci - r, C 2 - r) 

the four vertices of the square 5(r). By definition, ip = f at every vertex Vi: then, in 
the triangle Ti (and the same holds true in the triangle T 2 ) we have that the two affine 


functions cp and -0 satisfy \p — ijj\ = \ f — i/j] at every vertex, thus (4.6) gives 


Qt\/^ s 

Ilv^-V’I|L“(-S(r)) < ||/-V'|U-(-S(r)) < ||/“V'||L“(-S(2r)) < l2rR5+^- \\M\\ < 


and this, together with (4.6), implies the first estimate in (4.1). 


Concerning the second one, let us call Mi the constant value of Dp in Ti, and notice 

W2)-p{Vi))-[^{V2)-^{Vi)) 


that by (4.7) we get 

|(Mi-M)(ei) 


2 r 


£ 

^ 10 ’ 


and the same estimate holds true for 
V4. As a consequence, we have that t 


(Ml — M) (e 2 ) simply by checking the vertices Vi and 
le constant value of \Dp — D'ip\ in Ti is less than e/S, 
and the same holds true of course also in T 2 . In other words, \\Dp — T*t/'IU°°(*s(r)) < £^/5. 
Since DiIj is constantly equal to M in the square iS(r), by the definition of Lebesgue square 
we get 

\\Df - DpllL^isir)) = f \Df{z) - Dpiz)\dz < [ \Df-M\+ f \M - Dp\ 

Js(r) Js{r) Js{r) 


)S{r) 

4 

< 36r^6 H— r'^e < r'^e , 

5 


where the last inequality is true up to possibly decrease the value of <5: this gives the 


second estimate in (4.1) 


Let us now suppose that det M > 0. As a consequence, the image of S{r) under ip is 
a non-degenerate parallelogram, and then the image of S{r) under p is the disjoint union 
of two non-degenerate triangles, as soon as 11-0 — p\\L°°(sir))/r is small enough, depending 


on M; moreover, also the validity of (4.2) is obvious as soon as ||'0 —/||L°°( 5 ( 2 r))/F is small 


enough, depending on M and on e. Since (4.6) is valid for every w G 5(2r), we get (4.2) 


up to further increase the value of R and decrease the value of 5, again depending only 
on M and on e; thus, the proof is concluded. □ 
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Remark 4.4. Notice that, in the last estimate of the above proof the final values ofl/R 
and of 6 behave more or less as e multiplied by min |M(t)|/|t|, and the latter number 
is strictly positive exactly when detM > 0. This clarifies the need of the assumption 
det M > 0 m order to get (4-2). We can come to the same conclusion also directly by 
considering the claim of (f.2): we cannot hope it to be valid for the case when detM = 0; 
indeed, it is true that f is as close as we wish to an affine function, but this affine 
function is degenerate, hence the image of a small sguare around c is close to a degenerate 
parallelogram, which is a segment (or even a point if M = 0). And of course, knowing that 
the four vertices of a small sguare are sent very close to the vertices of a parallelogram does 
not even imply that the piecewise affine function (ps{x,r) is injective, if this parallelogram 
is degenerate! 


Remark 4.5. A guick look at the proof of the above lemma ensures that the constant 
(5(M, e) actually depends only on e, \\M\\, and det M: more precisely, S{M, e) is the mini¬ 
mum between a constant which continuously depend on e, \\M\\, and det M for det M > 0 
(found in the first part of the proof), and another constant which also depends continu¬ 
ously on £, \\M\\, and detM, but in the range detM > 0 (found at the end of the proof): 
this second constant explodes when detM —)■ 0. As a conseguence, 6 is bounded if e is 
bounded from below, ||M|| from above, and if det M is either d or it is bounded both from 
above and below (with a strictly positive constant). 


The crucial importance of the above general lemma comes from the fact that we can 
always apply it for small squares “almost centered” at Lebesgue points x. 


Lemma 4.6. Let 6 > 0 be given and let x be a Lebesgue point for Df. Then, there 
exists f = f(x,S) such that, for any r < f and any x G S{x,r/2), the sguare S{x,r) is a 
Lebesgue sguare with matrix M = Df{x) and constant 6. 


Proof. Let x E fl, M E and 5 > 0 be as in the claim. Since a; is a Lebesgue point, 
there exists f = f{x, 5) such that, for any r < f, one has B{x, 5r) C and 


-f Df{z) - M 

JB(x,5r) 


, (5 

dz<- 


(4.8) 


Let now x E S{x,r/2). We have S{x,3r) C B{x,5r) C Ll, and moreover (4.8) gives 
f \Dnz)-M\dz^fi^ f \Df-M\<,fiJ \Df-M\ 

Js{x,3r) Sbr^ Js(x.3r) JB(x.5r) 


36r2 JB{x,br) 

oO JBix.br) 


36r2 Js{x,3r) 
257r 


and by Dehnition 4.1 this means that S(x,r) is a Lebesgue square with matrix M and 
constant 6 . Thus, the proof is concluded. □ 


4.2. How to “move the vertices” of a grid. In this section we describe how to “move 
the vertices” of a grid in order to be able to control the average of \Df\ inside a square 
with the average of \Df\ on the boundary of the same square. To do so, we hrst introduce 
the following notation. 
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Definition 4.7. We say that the domain is an r-set if it is the finite union of essentially 
disjoint squares, all having side 2r and sides parallel to the coordinate axes. Any side of 
one of these squares will be called side of type A if both the endpoints are in the interior 
of side of type B if at least one endpoint is in dQ, but the interior of the side is inside 
n, and side of type C if the whole side is in dfl. Any vertex of one of the squares will be 
then called vertex of type A if it belongs to the interior offl, vertex of type B if it belongs 
to dVt but it is endpoint of at least one side of type B, and vertex of type C otherwise. 


For any small constant £ > 0, we will define a short segment or cnrve aronnd each 
vertex of type A and B. We start with the vertices of type A. 


Definition 4.8. Let be an r-set, £ <C 1 and let V = {Vi, V 2 ) be a vertex of type A. We 
call IfiV) the segment of length 2\/2er centered at V and with direction tt/A, namely 

Ie{V) = {(f"i +t,V 2 + t)-. \t\ < er} . (4.9) 

Notice that all the segments hiV) lie inside and they do not intersect with each other. 


The main resnlt that we prove ensnres that the average of \Df\ in a side of a sqnare 
can be always estimated with the average of \Df\ in the whole sqnare, np to move the 
two vertices in the corresponding segments fi. Actnally, in order to be able also to treat 
Lebesgne sqnares with det M = 0, we will in fact estimate \Df — M\ instead of \Df\ for 
some matrix M: then, we will apply this resnlt with M = 0 for all the non-Lebesgne 
sqnares, and with M = Df{x) for Lebesgne sqnares “almost centered” at a Lebesgne 
point X. Unfortnnately (bnt the reason is qnite evident) the estimate mnst explode as 
1/e; however, this will not be a problem for onr constrnction. 


Lemma 4.9. Let VL be an r-set, let AB be a side of type A, and let M G be a matrix. 
Calling TZ (2 Q the union of the 6 squares of the grid having either A or B (or both) as 
one vertex and 

T{A,B,M) = {{x,y)eIfiA)xUB)-. f \Df-M\dn^>-f |D/- M| , 

Jxy er Jn > 

one has 


HiQxGA(A): H^{[yeIfiB)-. (a;,i/)er}) 


> 


wiMmrw _ wiMA)) 


7 


< 


(4.10) 


Proof. Let ns snppose, jnst to fix the ideas, that the side AB is horizontal, as in Fignre[TTj 
snppose also, for the moment, that M = 0. Let then 77o C 77 be the small parallelogram 
-dark shaded in the fignre, while 77 is light shaded- whose fonr vertices are the endpoints 
of the segments A (A) and IfiB). A simple change of variable argnment, together with 
the fact that 77o C 77, ensnres that 

[ [ [ iDfldPL^dydx^Ser [ \Df{z)\d'H‘^{z)<8er [ \Df{z)\d'H'^{z). 

Jxeh{A) Jy&Ie{B) Jxy JUo Ju 

On the other hand, writing F = F(A, B, 0) for shortness, we also have 


25 
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A 

TZo 






B 


Figure 11. The rectangle TZ, the side AB, and the two segments Ie{A) and 


hence we dednce 




from which (4.10) immediately follows. 


To show the general case with M 7 ^ 0, it is enongh to apply the above argnment to the 
fnnction /(x) = /(x) — Mx: of conrse th e set r( 74 , B, M) coincides with the set r(y4, 0) 
corresponding to the fnnction /, hence (4.10) follows also in the general case. □ 


The above resnlt will be nsefnl in order to treat the internal sqnares, bnt we have to 
take care also of the bonndary sqnares. Let us then extend the dehnition of the segment 
hiy) to the vertices of type B. 

Definition 4.10. Let Vt he an r-set and let V G dQ be a vertex of type B. IfV is vertex of 
exactly two squares of the decomposition, and these two squares are adjacent (as for Vi in 


Figure 12, left), then we call leiV) the segment of length 2er on d^l centered at V. IfV 
is vertex of three squares of the decomposition (as for V 2 in the figure), then we call leiV) 
the union of the two segments contained in dQ, of length er, having V as one endpoint. 


-H 

H- 

-F2 





“1-- 

F3 










S2 

□ 

^3 


Figure 12. Left: some squares and vertices Vi near the boundary of an 
r-set n, and the corresponding IsiVi). Right: dehnition of 7y and Ty. 


Notice that we have dehned Ie{V) only for the 
instance not for points as V 3 of V 4 in Figure [T^ We 


Lemma 4^ for sides of type B. To do so near points 
last simple dehnition. 


vertices of type A and B, thus for 
want now to extend the validity of 
like V 2 in Figure [T^ left, we need a 


Definition 4.11. Let 11 be an r-set and let V G df2 be a vertex of three squares of the 
decomposition, say iSi, S 2 , S 3 . We call Tv the right triangle having right angle at V, two 
sides of length r/2, one horizontal and one vertical, and being not contained in 11, and for 
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i = 1, 2, 3 we call S~ the square contained in Si, having one vertex at V, and side r 12. 
Then, we let he the obvious piecewise affine homeomorphism between U5^UiS3" and 
U 53 " U7y, which is hi-Lipschitz with constant 2. Finally, we call fi"*" the union of 
with all the triangles 7y for vertices V as above, and we call —)■ the piecewise 

affine homeomorphism which coincides with around every vertex V, and which is the 


identity outside (see Figure 12. right). Notice that ^ is the identity in a r-neighborhood 


of all the internal squares of the decomposition, and it is globally 2-biLipschitz. Finally, 
for every x, y E fl such that the segment ^(a;)^(?/) is contained in we call xy the 
counterimage, under 'P, of this segment (which is of course a piecewise linear path). 


We can finally generalize Lemma 4.9| for all the sides of type B; it will be enough to 
limit ourselves to consider the simpler case M = 0 . 


Lemma 4.12. Let VL be an r-set and let AB F Q be a side of type B. Calling TZ F Q the 
union of the squares of the grid having either A or B (or both) as one vertex and defining 

100 

Jn 

OTIC has 

(j:,i/)er})>— 


ry,B) = {(i.!/) 6 CO)xC(B): L\Df\dW>—l \Df\dHA. 

Jxy er Jn > 


Proof. This is a very simple generalization of Lemma 4^, we just have to consider a few 
possible cases, all depicted in Figure [T^ Without loss of generality, we can assume that 
B G dQ and that the segment AB is horizontal. If B belongs to three squares of the 
decomposition, then there are three possible subcases. First of all, A can be inside 12 
(this is the hrst case depicted in the hgure); second, A can also belong to three squares 
of the decomposition (depending on how these squares are, this is the second or the third 
case in the hgure); last, A can belong to two squares, which are then necessarily the two 
squares having AB as a side (this is the fourth case in the hgure). Otherwise, B can 
belong to exactly two adjacent squares of the decomposition, and then again A can either 
be inside 12 (hfth case in the hgure) or in the boundary of 12 : in this latter case (the sixth 
and last one in the hgure) A must also belong only to the same two squares to which 
B belong, since otherwise we fall back into an already considered case. To clarify the 


situation. Figure 1^ does not show the situation in 12, but directly in 12+. 

The proof is now almost identical to the proof already done in Lemma [4.10[ Let us 
call again TZq the quadrilateral having as vertices the endpoints of ^/(/^(A)) and ^/(/^(i?)); 
this quadrilateral, depicted in the hgure for all the possible cases, belongs to 12+. Notice 
that, depending on the case, IZ can be done by 2, or 3, or 4, or 5 squares, and the hgure 
always shows only (the image under T of) these squares. To conclude the proof we only 
have to keep in mind that we are interested in what happens in the real domain 12 , not 
in the simplihed domain 12+. However, we can use the map T to move the situation from 


12 to 12+; then, we notice that the very simple calculation done in Lemma 4.10 works 
perfectly in the new situation; and hnally, we use to get back to the case of 12. The 
only detail which changes, since T is 2-biLipschitz, is that the constant 25 in the old 
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Figure 13. The six possibilities in Lemma 4.11 


definition of T for internal sides becomes 100 for the new definition of T for sides touching 
the boundary. □ 

Let us now fix a matrix M = M{A,B) for any side AB of type A, and write for 
brevity r(y4, B) = r(y4, B, M{A, B)). Thanks to the above results, we have defined a set 
of “bad pairs” {x,y) G T{A,B), where “bad” means that in the segment xy (or in the 
curve xy) too much derivative is concentrated. We can now find a selection of points 
Xy G leiy) for any vertex V, so that for any side AB the pair {xa_,xb) does not belong 
to T{A,B). 

Lemma 4.13. Let hi be an r-set, and let us fix a matrix M{A, B) for any side of type A. 
It is possible to select a point xy G hiV) for any vertex V of type A or B in such a way 
that, for every side AB of type A or B, one has {xa,xb) ^ r(A,il). 


Proof. We will argue recursively. First of all, let us enumerate all the vertices of type A 
or B as Vi, V 2 , ... Vat, for some A^ G N. Then, we aim to show that it is possible to select 
by recursion points x* = xy^ in every A (14) in such a way that, whenever 1414n is a side 
of type A or B, the point Xi is chosen in such a way that 




yelfiVj: (x„|/)Gr(V),VA}) < Am>t, 


(4.12) 

(.Tj, X^n) 

Notice that, since we will define the points recursively, then the above requests make 
sense: in particular, if m < z then the point Xm has been already chosen when we have 
to choose Xi- Of course, if we can select all the points Xi according to (4.12), then we are 
done: the thesis will be simply given by the second property in (4.12), but the first one 
is essential to let the recursion work. 

Let then 1 < i < A^, and let us suppose that the points Xj for j < i have been already 
chosen according to (4.12); let n~ (resp., n'^) be the number of the indices j < i (resp., 
j > i) such that is a side of type A or B. By (4.12) applied to the indices j < i, we 
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know that the points x G leiVi) such that {x, Xj) G r(i^, Vj) for some j < i corresponding 


to a side ViVj cover a portion at most n /5 of Is{Vi). On the other hand, by Lemma 4.12 
the points x G leiVi) such that 


nH\yeh{V,)-. {x,,y)eT{V,,V, 


> 




for some j > i for which VjVj is a side of type A or B cover a portion at most n'^/5 of 
leiVi). Since of course n~ + n'^ < 4, we can pick a point Xi G A (hi) for which none of 
the above problems occur, hence by definition this choice fulfills (4.12). The recursion 
argument is then proved, and the proof is concluded. □ 

The last goal of this section is to define an approximating function / of / on the grid 
given by the boundaries of the squares. First of all, let us give the definition of “grid” 
and “modified grid”. 

Definition 4.14. Let hi be an r-set, and for any vertex V of type B let V he a given point 
xv in leiy); instead, let V' = V for any vertex V of type A or C. We call grid the union 
Q of all the sides AB of the squares of the decomposition, while the modified grid is the 
union Q of all the “modified sides”, that is, the piecewise linear curves A'13'. Notice that . 


if AB C dVL, it might happen that the curve A'B' has not been defined in Definition f.ll 


if this is the case, we simply denote by A'B' the shortest curve in hi connecting A' and B': 
notice that this shortest curve lies entirely inside dQ, and that actually this minimizing 
property f or A'B ' is true also for the sides AB C dfl where A'B' was already defined in 
Definition f.ll. For every square S of the grid, we call S the union of its modified sides. 


Observe that the grid Q coincides with the grid Q, except near the boundary of hi; 
analogously, the piecewise linear curve A'B' is nothing else than the segment AB, if it is 
a side of type A. Notice that both the grid and the modihed grid contain the boundary 
of hi. We give now our last dehnition of a map / on Q. 

Definition 4.15. Let hi be an r-set, and for any side AB of type A fix a matrix M = 
M(A, B ). Le t the points xy G leiV) for the different vertices V of type A or B be as in 
Lemma f.lS . We define the function g : Q ^ as follows. For any side AB of type A 
or B, we define g on the curve A'13' as the reparameterization, at constant speed, of the 
function f on the curve xaXb; moreover, we let g = f on dQ C Q. 

In the above definition, it is important not to confuse the points xy with the points 
V'\ according with Dehnition 4.14 V' = xy if D is a vertex of type B, while V' = V \iV 


is a vertex of type A or C. In particular, if AB is a side of type A, then A'B' is simply 
the segment AB, hence g on the segment AB is the reparameterized copy of / on the 
segment xaXb = xaXb- We conclude this section with an estimate for the function g. 

Lemma 4.16. Let hi be an r-set, and let the matrices M = M{A,B), the points xy and 


the function g : Q ^ be as in Definition f.lS. Then, for any side AB of type A 


calling u the unit vector with direction AB and IZ again the union of the squares of the 










38 


STANISLAV HENCL AND ALDO PRATELLI 


grid having either A or B as one vertex, we have 

25 


[ \Dg(t)-M-u\dt< — f \Df - M\d'H‘^+ ll\\M\\er. (4.13) 

JAB er Jn 


Instead, for any side AB of type B, we have 

j^\Dg{t)\dt< — 
JA'B' er 


'n 


\Df\dn\ 


(4.14) 


Proof. Let us start with a side AB of type A, and call for brevity x = xa and y = xb- 
By Lemma 4.9, we already know that 

[ \Df(s) - M\ds <— f iDf-MldH^ (4.15) 

Jxy ST J'R, 

Recall now that, by definition, the function g on the segment AB is simply the reparam¬ 
eterization of the function / on the segment xy. Define then A = AB/xy, and call v the 
unit vector with direction xy. notice that by construction 


l-2e<A<l + 3e, 

As a consequence, by a change of variable we obtain 

Df{s) ■ h 


<2e. 


'AB 


\Dg{t) — M ■ i>\dt = \ 


/xy 


A 


-M-u 


ds = 


/xy 


Df{s) ■ u — M ■ Xu 


ds 


< f {Df{s) — M) ■ u ds + f \\M\\ \Xu — u\ ds 

-I xy J xy 

< [ \Df{s) - M\ds + 5\\M\\exy, 

J XV 


thus recalling (4.15) we get (4.13). 

Let now AB be a side of type B, and call again for brevity x = xa and y = ys- In 
this case, by Lemma 4.12 we already know that 

/* 1 DO r 

l\Df\dn^<— \Df\dn\ (4.16) 

J xy ET JtZ 

Now, by definition of g we have that the image of A'B' under g coincides with the image 
of ah/ under /, hence the lengths of the two curves coincide, which means 


i\Df\d'H^= A |D(?(f)| 

Jxy JA'B' 


dt . 


Hence, (4.14) directly follows from (4.16), and the proof is concluded. 


□ 


4.3. How to modify / on a grid. In this section, we show how to modify a function 
on a one-dimensional grid; more precisely, we take a generic function defined on a grid, 
and we modify it in order to become piecewise linear. We have to do so because both 
our big results, namely Theorems 2T and 3T, need a function which is piecewise linear 
on the boundary of a square. We start with a rather simple modification, which we will 
eventually apply to the “bad” squares and to the “good” squares corresponding to a 
matrix with det M = 0; this construction is reminiscent to the one made in imi. where 
the situation was more complicated because also the inverse should be approximated. 
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Propositio n 4.17 . Let be an r-set, let Q and Q be a grid and a modified grid according 
to Definition \4-. and let g : Q -f-be a continuous, injective function which is piecewise 
linear on dQ. Then, there exists a piecewise linear and injective function g : Q ^ such 
that g = g on dfl, giV') = giV') for every vertex V of the grid, and for every side AB of 
type A and every matrix M one has 


[ \Dg{t) — M ■ ul dt < [ \Dg{t) — M ■ ul dt, 

JAB JAB 

while for every side AB of type B one has 


(4.17) 


/- \D!l{t)\dt < /- \Dg(t)\dt. 

Ia'B' Ja'B' 


(4.18) 


Moreover, on each curve A'B' the function g is an interpolation of finitely many points 
of the curve g{A'B'). 


Proof. We define the map g in two steps, namely, we work first aronnd the vertices, and 


then in the interior of the sides. Fignre depicts how the constrnction works. 

Step I. Definition of g around the vertices. 

Let ns start by selecting a vertex V of type A, that is, V belongs to the interior of hi. There 
are then fonr sides of the grid for which V is an endpoint, and we can call Vi, 1 < ^ < 4 
the four other endpoints of these sides. Since g is continuous and injective, the quantity 


inf \g{x)g{V) : t G ^ \ (J W 


2=1 


is strictly positive: let us then select a small radius p = piV) > 0, much smaller than 
this quantity. Hence, by definition, the ball B{g{V),p) intersects the image of the grid 
Q under g only in points of the form g{x) for x belonging to one of the four sides VVi. 
On the other hand, for each z = 1, ... , 4 there is at least a point x G VVi such that 
g{x) belongs to dB{g(y), p). Let us then call Vfi the last point x of the segment VVi for 
which g{x) G dB{g{V), p), where “last” means “the farthest from V”. We define then the 
function g, on each of the four segments VVA, simply as the linear function connecting 
g{V) with 

Let us now consider a vertex V of type B, that is, V belongs to dfl but there is some 
side of the grid, contained in the interior of O, of which V is an endpoint: let us call j 
the number of such sides, and notice that by construction j is either one or two (keep in 
mind Figure [I^. We argue then similarly as before: we call Vj, for 1 < i < j the other 
endpoints of these internal sides, and we consider the strictly positive quantity 

inf [g{x)g{V) : x e g \ A dfl)] . 


This time, we will define p = piV) not only much smaller than the above quantity, but 
also so small that g{dfl) fl B{g(y'),p) is the union of two segments (this is surely true 
as soon as p is small enough, since g is piecewise linear on dfl). Exactly as before, for 
1 < * < j we call Vfi' the last point x G V'Vj such that g{x) G dB{g(y'), p): up to further 
decrease piV), we can also assume that the portion of the piecewise linear curve V'Vf 
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connecting V and is simply a segment. Then, as before we define g on each of the j 
segments as the linear function connecting giV') with giV^). 



Figure 14. The construction in Proposition 4.17 the points A, B and 
are in while D, E, F and G are inside the image of g inside 
thicker. 


C 

is 


Step II. Definition of g inside the sides. 

Up to now, we have dehned ^ on a neighborhood of each vertex of type A and B; moreover, 
g is already automatically dehned on dfl, since it must he g = g on dQ. Therefore, 
to conclude we have to dehne g on the remaining portion of Q. By construction, this 
remaining portion is a hnite and disjoint union of internal parts of sides of type A or B; 
more precisely, for every side AB of type A there is a segment A^B~ CC AB where g 
has to be dehned, while for every side AB of type B g has still to be dehned in some 
piecewise linear curve A+B~ CC A'B'. 

Let us hrst consider the case of a side AB of type A: the function g has been already 
dehned in the segment AA'^ (resp., B~B) as the linear function connecting g{A) and g{A'^) 
(resp., g{B~) and g{B)), and moreover the points g{A~^) = g{A~^) and g{B~) = g{B~) are 
in the boundary of the disks B{g{A),r{A)) and B{g{B),r{B)) respectively. We have to 
dehne g in the segment A'^B~ , and this must be a piecewise linear curve connecting g{A~^) 
and g{B~). Observe that g, in the segment is already a curve connecting g{A~^) 

and g{B~), the only problem being that it is not necessarily piecewise linear. We can then 
select many points Pq = A+, Pi, P 2 ,... , Pn = B~ in the segment A+P“, and dehne g 
in A'^B~ as the piecewise affine interpolation of these values (that is, g{Pi) = g{Pi) and g 
is linear on each PjPj+i). A very simple geometric argument (similar to [101 Lemma 5.5], 
but much easier) shows that, by choosing carefully many points, the map g in A'^B~ 
is injective, it never crosses the two disks B{g{A), p{A)) and B{g{B), p{B)), and its L°° 
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distance with g is much smaller than 


inf [g{x)g{y) : x G A^B , y eQ\ AB^ . 

Let us now consider a side AB of type B; in this case, g is piecewise linear in the two 
segments A'A'^ and B~B', and we have to extend the dehnition to the piecewise linear 
curve A+B~. This can be done exactly as we just did for a side of type A, the only 
difference is that, by doing the interpolation, the points Pi must include all the extremes 
of the segments forming the curve Apart from that, nothing else changes, thus 

the dehnition of g is hnally concluded. 

Step III. The properties of g. 

To conclude the proof, we just need to check that g fulhlls all the required properties. 
The fact that g = g on dQ and at every vertex is true by construction, as well as the fact 
that g is an int erpola tion of hni tely many points of the curve g{A'B') on every side A'B' 
of Q. To check (4.17) and (4.18), we just have to keep in mind that /as \Dg\ is the length 


of the curve g on the segment AB, while Iab \Dg\ is the length of the (piecewise affine) 
curve g on AB. But any interpolation of points of a curve is shorter than the curve itself. 


thus (4.17) follows immediately for the case M = 0, and the very same argument with 


A'B' in place of AB ensures also (4.18). 


To show (4.17) when M 7 ^ 0, let CD C AB be a segment where g is linear and satishes 


g[C) = g{C), g{D) = g{D): then 

rD 

I \Dg{t) — M ■ z/| dt = 


ic 


rD I'D 

/ Dg{t)—M-vdt = / Dg{t)—M-vdt 

Jc Jc 


< 


[ \Dg(t) — M ■ u\ dt, 
Jc 


thus adding over the segments where g is linear we get (4.17) and the proof is hnished. □ 


We conclude this section with the following generalization of the above result, rather 
technical but very useful to obtain Theorem 1.1[ and whose proof is actually nothing else 
than a straightforward modihcation of the previous one. 

Proposition 4.18. Let C = UfL^AiBi be a finite union of closed segments in and 
let Co = with Nq < N, he a selection of some of them. Let : C —)■ &e 

a continuous one-to-one function, and let rj > 0 be given. Then, there exists another 
continuous and one-to-one function g -. C such that g = g at every endpoint of each 

segment, {g 7 ^ g} is contained in the rj-neighborhood of Co, g is piecewise linear on Co, 


and for every side of C the estimates (4-1I) cind (f.lS) hold. 


Proof. The proof can be done almost exactly as in Proposition |4.17[ First of all, up to a 
subdivision of some of the segments, we can assume that every two segments are either 
disjoint, or they meet at a common endpoint. Then, we start dehning ^ on all the 
sides both whose endpoints are not contained in Co- Further, we consider any of the 
remaining vertices, say V. We can select a small p such that the ball B{g(y),p) contains 
only points of the form g{x) for x belonging to one of the (hnitely many) segments VVi 
having V as an endpoint; up to decrease p, we can also ensure that \x — V\ < p for any 
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such X. We define then the points exactly as in Step I of Proposition 4.17, and we 


let g be linear on each segment VVj^: the continuity and injectivity of g up to now is 
then clear. In the portions of the segments where g has not yet been defined, we can 
then define it in two different ways: inside the segments which belong to Cq, we define a 


piecewise linear g exactly as in Step II of Proposition |4. 1 7; inside the other segments, we 
simply let g = g. 

By construction and arguing as in Step III of Proposition |4. 17 we can then immedi¬ 
ately observe that the function g is as required. □ 


4.4. The proof of Theorem |1.1[ We are finally in position to give the proof of Theo¬ 
rem 1.1, which will come by putting together all the results that we got up to now. For 


the reader’s convenience, we split the proof in some parts. The first one is a very peculiar 
case, namely, when is an r-set and the function / is already piecewise linear on the 
boundary; nevertheless, most of the difficulties are contained in this part. 

Proposition 4.19. Under the assumption of Theorem [O assume that in addition Ut is 
an r-set, and that f is continuous up to dQ and piecewise linear there. Then, for every 
£ > 0 there exists a finitely piecewise affine homeomorphism such that 


life — /llwm + ll/e — /||l“ < £, 


fe = f on dVt. 


(4.19) 


Proof. The idea of the construction is rather simple: we divide the squares in four groups, 
namely, the Lebesgue squares with positive determinant, the Lebesgue squares with M ^ 
0 but zero determinant, the Lebesgue squares with M = 0, and the other ones. Inside the 


first squares we can substitute / with (ps{c,r) and rely on Lemma 4.3, for the second ones 


we will use Theorem |3.1[ and for the third and fourth ones Theorem |2.1[ However, to treat 
each square separately, we need to take care of the values on the boundaries of the squares: 


on one hand, they must be piecewise linear, in order to allow us to apply Theorems 3.1 


and 2.1, and this will be obtained thanks to Proposition 4.17 but on the other hand, any 
two adjacent squares must have the same boundary values on the common side, and this 
will require same care. Let us then start with the proof, dividing it in several steps. 

Step I. Definition of the constants Si and of the sets Ai, A 2 and A 3 . 

First of all, we have to take five small constant et for 1 < i < 5. More precisely, Si is 
a small geometric constant (for instance, ei = 1/10 is enough); instead, the constants 
£5 £4 £3 <C £2 ^ £ will depend on the data, that is, fl, / and e. More precisely, 

since / G IF^’^(f 2 ), we can select £2 ^ 1 so small that 


/l^/l< — 

I a' - 5477 


VH C : |H| < £ 2 , 


(4.20) 


where 77 is a purely geometric constant, which we will make explicit during the proof. 
Then, we define £3 -C £2 in such a way that 


X G 72 : 0 < \Df{x)\ < 83 or \Df{x)\ > — or 0 < det(71/(x)) < 83 

^3 


<|. ( 4 . 21 ) 
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This estimate is surely true as soon as £3 is small enough, depending on £ 2 , op ^5 and on 
/. Now, we define the following two subsets of the matrices M G 


< ||M|| < — , det M > eaj , < \\M\\ < — , det M = o| , 

^ (S3 * ^ Sq 


1 

^3 


which of course depend only on £ 3 . Finally, we let £5 -C £4 £3 be so that 


^5 ^4 ^5 ^ £ 

£4 £3 ~ GK\il\ ’ 


£5 <c £4^3 , 


£4 < ^3, 


(4.22) 


and we set <5 = 6 ( 83 , £ 5 ) as 


S = min {S{M, 85 ), M e M+UM^} , (4.23) 

Observe that S really depends only on £3 and £5 


4.3 


4.5 The last constant to select is r: indeed, O is 


where 5 are the constants of Lemma 
by construction, as observed in Remark 
an r-set, but then we can regard it as an r/H-set for every FT G N; as a consequence, we 
can now change the value of r, making it as small as we need: in particular, we let r be 
so small that 

{a: G O : f(x,(5) < r}| < ^ , ^ 


rP{n) 


rP{f{Q)) < 


(4.24) 


4.6 


66 K ’ 

for every a; G O which is a 
where P{A) is as usual the 


where the constants f{x, 6 ) have been defined in Lemma 
Lebesgue points for Df (so, for almost every point of O), anc 
perimeter of A, that is, P}{dA). 

Having fixed all the constants £*, and having also chosen the final value of r, we can 
now enumerate the squares of the grid as 1 < i < N, and we subdivide these squares 
in four groups, namely. 


Ai = CC n : is a Lebesgue square with matrix Mi E Ai~^ and constant , 

A 2 = CC hi : is a Lebesgue square with matrix M* G and constant (jj , 

*^3 = |‘5j CC is a Lebesgue square with matrix Mi = 0 and constant , 

v 4.4 = liSj : Si ^ Ai U *4-2 U 4.3 


We aim now to show that most of the squares belong to the first three groups. In fact, 
let us consider the total area of the squares containd in A 4 . The union of those which 
touch the boundary of hi has of course an area smaller than rP{Q). Let us instead 
consider a square S{c, r) G Ai compactly cont aine d in hi: this means that, for every point 
X G S{c,r/2), either we cannot apply Lemma 


4.6 


with constant <5 to x (thus, r(x, 5) < r), 
or X belongs to the set in (4.21). Since the area of S{c,r) is four times greater than the 


area of S{c,r/2), by ( |4.21[ ) and (4.24) we deduce 

1^ *^ 4 } < rP(f2) +4^ {x G 12 : f(x, 5) < r] + 




(4.25) 


Step II. Squares in Ai and A 2 never meet. 

Let us now show that squares in Ai and A 2 never meet, that is, no vertex of a square in 
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can be vertex also of a square in A 2 '- this will come as a simple consequence of the 
L°° estimate (4.1) in Lemma 4.3 Indeed, assume for simplicity of notations that iSi = 
5(ci,r) and S 2 = S(c 2 ,r) have some common vertex, and that 5i G Ai^ S 2 G A 2 - Then, 
Lemma [4^ provides us with two affine functions " 01 , "02 satisfying D'lpi = Mi, = ^ 2 , 
each Si being a Lebesgue square with matrix Mi. Let us now call 


‘^ 3/2 S 


Cl + C2 


) C 5(ci,2r) n 5 (c 2 , 2 r) 


so that by (4.1) and recalling (4.23), we have 


IIV'l -'02 ||l°°( 53/2) < 11^1 “ /|U°°(5(ci,2r)) + ||^2 “ /||l,“(5(c2,2r)) < 


re^ 


Since det M 2 = 0, by construction we can hnd two points x, y E S ^/2 such that \y — x\ = r 
and 'ip 2 {y — t) = 0. The last inequality, keeping in mind the dehnition of ^ 1 , yields then 

re^ det Mi , 

^ > m[y -x)\> \y-x\> . 


Since this is in contradiction with (|4.22|), we have concluded the proof of this step. For 

(4.26) 


later use we underline that, more in general, we have proved what follows: 

C5 


V Sa, Sb G *4-1 U A 2 adjacent, one has \\Ma — Mb\\ < — . 

5 

Step III. A temptative modified grid Qi and of a temptative function gi : Qi ^ . 

In this step, we dehne a modihed grid Qi and a function gi on it. To do so, we simply 
have to choose a matrix M{A,B) for every side AB of type A of the grid Q] once done 


this, we get hrst the points xy G hfiV) from Lemma 4.13 applied with £4 in place of e, 
then the modihed grid Qi from Dehnition 4.14 and hnally the function gi '. Gi ^ from 
Dehnition 14.151 

The matrices M{A,B) will be dehned as follows. For any side AB of type A, let 
for a moment Sa and Sb be the two squares of the grid having AB as a side: then, if 
neither Sa nor Sb belong to A 2 , we let M{A, B) = 0; if Sa G A 2 but Sb 4- *^ 2 , then we let 
M{A, B) = Ma, and analogously if Sa ^ A 2 and Sb G A 2 we let M{A, B) = Mb, hnally, if 
both Sa and Sb belong to A 2 , then we set M{A, B) to be arbitrarily one between Ma and 
Mb', it makes no diherence which one we choose, since by (4.26) we know that M^ Mb. 

We want now to evaluate jQ^\Dgi\ for some of the modihed squares S (recall Def¬ 
inition 4.14). Let us start by taking a square S = S(c,r) G 4,3 U A 4 , and let us call 


iS’*' = S (c, 3r) n D the union of the nine squares around it (more precisely, of those which 
are inside D). Take any side AB C dS of type A or B, and observe that the union TZab 
of the squares touching A or i? is contained in If AB is of type A but M{A, B) = 0, 
or if AB is of type B, we can apply Lemma 4.16 (using ( |4.13 ) or (4.14) if AB is of type 
A or B respectively) and get 
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Instead, if AB C dS ft clff, then of course 


\Dgi\dH^ = f \Df\ dU^ = f \Df\ dV} . 

J A'B' JXjUXg J X iXnndQ 


Xj^Xgndfl 


Since the boundary of S is nothing else than the union of its four modihed sides A'B', 
adding the last two estimates for the four sides of dS we get 


[^\Dg^\dn^< — 
Ids e^r . 


/ \Df\dn^+ iDfldn^ WSeA- 

's+ Jds+nan 


■3,4 ) 


(4.27) 


being 

*^ 3,4 = 1*5 G ^3 U ^4 : M (A, 5) = 0 for each side AB of S [ 

Let us now consider a square S = S(c,r) G A 2 , and notice that by dehnition it is 
compactly contained in so 5 = 5 and A'B' = AB for any of its sides. Let AB be one 
of those sides, and notice that \\M{A^B) — M\\ < £ 5 / 5 , being S a Lebesgue square with 
matrix M and constant 5\ indeed if the other square having AB as a side is not in A 2 , 


then M{A,B) = M, and otherwise the inequality is given by (4.26). As a consequence, 
again (4.13) of Lemma 4.16 and the dehnition of give 

2 


f \Dgi-M-lyldH^ < [ {Dg^ - M(A, B) ■ u\dn^ + ^ re 5 

JAB JAB 5 


< 


25 

e:4r 

25 


/ \Df - M{A,B)\dn^ + ll\\M\\e^r + - re^ 

In 5 

< —4 [ \Df — M\ d'H'^ + 120r — + Hr — H— re^ 

S^T JlZ £4 £3 5 

< —^ / \Df — Ml dl-L^ + 120r — + Hr — H— re^ 

e^r JS+ £4 £3 5 

< 1020r — + Hr — + - re^ , 

£4 £3 5 


(4.28) 


4.1 


together with the fact that S < 


where in the last inequality we have used Dehnition 
5{M, £ 5 ) < £ 5 . ^ ^ 

Step IV. The “correct” modified grid Q and a second temptative function g 2 : G ^ . 

In this step, we dehne a second modihed grid and a second temptative function; the idea 
is to repeat almost exactly the procedure of Step III, but using the neighborhoods hfiV) 
instead of those fifiV). In fact, the presence of £4 is perfect for the squares in A 2 , since 


in (4.28) we only have small terms like 64/83 or 63 / 84 , instead, for the squares in A 3 UA 4 , 


the constant 64 in (4.27) is too small, and we would need something much larger than 
62 - Since there is no constant which is at the same time much larger than 62 and much 
smaller than £ 3 , we are forced to repeat the procedure. 

This time, let us dehne the matrices M'{A, B) = 0 for all the sides AB of type A, 
and let us consider a slightly modihed version of the intervals More precisely, we 

let 1'^.^iy) = /ei(H) for all the vertices V which are not in the boundary of some square of 
Ai or H 2 . Instead, if a vertex V = (Hi, V 2 ) belongs to the boundary of at least a square 
in Ai or A 2 (these two things cannot happen simultaneously, thanks to Step II), let us 
dehne IfiiV) as a translation of leiiV) of (± 2 £ir, ± 2 £ir), where the two choices of the 
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sign ± are done in snch a way that the whole interval is inside a sqnare of Ai or A 2 ] for 
instance, if V is the lower-left corner of a sqnare in Ai (or ^ 2 ), we can set 


4 (^) = {(^1 + t, V2 + t): Sir <t< 3eir} , 


compare with (4.9). If V is corner of more than one sqnare in Ai or A 2 , then we let 
be inside one of them arbitrarily, this will not make any difference. Fignre Ifo shows an 
example of a portion of an r-set where eight intervals are depicted and the 

shaded squares are those in Ai (or ^ 2 )- Notice that the intervals and could 

be inside each of the two shaded squares, in this example we have put the hrst interval 
inside the above square and the second interval in the below one. 



Figure 15. The intervals in Step IV. 


After a quick look at the proof of Lemma |4.9| , it is evident that it works perfectly 
even with the intervals in place of the /^^(V): indeed, in that simple proof we 

just used that the internal intervals are all of length 2y/2er, with direction at 45°, and 
placed very close to the vertices. As a consequence, we obtain the points Xy G 
from Lemma 4.13, the modihed grid Q from Dehnition 4.14, and the function g 2 Q ^ 


from Dehnition 4.15 The modihed grid that we obtained now is the “correct” one, and 
we will use the function gi (resp., ( 72 ) around squares in A 2 (resp., ^3 and A 4 ). The very 
same calculations of last step work also for this new case, just substituting the constant £4 
with El- In particular, since this time M'{A,B) = 0 for all the sides, the estimate (4.27) 
is true for every square in ^3 U .4,4, so we can rewrite it (with Ei in place of £4) as 


j^\Dg2\dV} < — 

hS EiT 


's+ 


\Df\d'H^ + 


' dS+r\dEl 


\Df\dV} ySeAsUAi 


(4.29) 


Step V. Definition of the function g^ : Q ^ . 

We are now in position to dehne a function g^ on the grid Q introduced in Step IV. This 
function will behave almost correctly on the whole grid, its only fault (which will be solved 
in next step) being not to be piecewise linear. As we already observed, we would like to set 
93 = Vs OLi the boundary of any square S & Ai, g^ = gi on the boundary of the squares 
in 4 - 2 , and 773 = g 2 on the boundary of squares in 4,3 or 44 ; of course, this is impossible 
because the function g^, would then not be continuous and injective. As a consequence. 
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we use the above overall strategy to define g^, but with some ad hoc modification where 
squares of different types meet, so to get continuity and injectivity. 

Let us start with the easy part of this definition. For every side AB which is in the 
boundary of some square of Ai, we define on AB as the linear interpolation which 
satisfies (73(A) = /(A) and gs^B) = f{B): as a consequence, (73 = ips on dS for every 
5 G Ai, where (ps is given by Definition | 4 . 2 [ Second, for every side AB which is in the 
boundary of some square in A2, we let 773 = 771 on AB: recall that vertices of squares in 
Ai and vertices of squares in A2 are disctinct by Step II. Finally, for every side AB such 
that neither A nor B are vertices of squares of Ai or of A 2 , we let 773 = 772 on A'B\ where 
the points A' and B' are those given by Step IV. Notice that, up to now, the function 773 is 
continuous and injective: this comes as a straightforward consequence of the L°° estimate 
around squares in Ai, and of the fact that 773 is a reparameterization of / on different 
segments around squares in A2 or A3 U A4. Nevertheless, 773 has still not been defined in 
the whole Q. 

Let us then consider a side AB such that 773 has not yet been defined on A'B': by 
construction, this means that at least one between A and B is vertex of a square in Ai 
or A2, thus in particular AB is not in the boundary of D, and both the squares of the 
grid having AB in the boundary belong to A3 U A4. We aim now to define g^ on A'B'. 
To do so, let us keep in mind that we would like to set gs = g2, and observe also that 772, 
on A'B', is nothing else than the reparametrization, at constant speed, of the image 70 of 
some piecewise linear curve under /. Our idea is to define g^, on A'B', again as the 

reparametrization at constant speed of some modification 7 of 70. In particular, 7 and 70 
will coincide in their big “internal” parts, the difference being only near the endpoints of 
these curves. 



Figlfre 16. The definition of 773 in Step V if 5 is in some square of A2. 


For simplicity, let us start assuming that A is not vertex of squares in Ai or A2, 
and B is vertex of at least a square in A2 (by construction and recalling Step II, B is 
then vertex of either one or two squares in A2, and of no square in Ai). By Step IV, 
we know that the interval I'^^{B) is entirely inside a square S of A2 which has B as a. 


vertex. As in Figure 16 (left), let us call C ^ B the vertex of S such that there is a 
square having both A and C as vertices. The function g^ has already been defined, on 
BC, as the reparametrization of the image of a segment xbXc under /; by definition and 
by construction, the segment xbXc and the curve x'j^x'^ meet in some point P near B] 
notice that B & Q, and then the curve x'j^x'^ is actually a segment, except possibly in a 
small neighborhood of We are now ready to define the curve 7: first, we define 7 as 
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the image, under /, of the union between the part of from xa to P and the segment 
Pxb- Then, since 7 and gzixsxc) have of course a part in common, we let 7 be a slight 
modihcation of 7 which intersects gsixsxc) only ai B = /{xb)- Of course, we need to 
modify 7 only between f{P) and B: this modihcation can be done as the enlargement of 


Figure 16 (right) shows, and it works exactly as in Step 10 of the proof of Theorem 2.1 


in particular, the length of 7 is as close as we wish to the length of 7 . By dehnition, we 
have 

^'(7) < [~\Df\ + f \Df\ = [ \Dg,\ + / \Dg,\ . (4.30) 

Jx'x' Jxb^c JA'B' JbC 


Thanks to (4.28), calling for a moment M the matrix associated to 5, we know that 


[ \Dgi\< [ \Dg^-M ■ u\dn^+ 2 r\M\< (1020-+ 11 -+ le5)r+- [ \Df\, 

JBC JBC \ £4 Ez o J r JS 


and then (4.30) becomes 


V-\l)< j^\Dg2\ + {lD2D^- + ll- + \eX + - j \Df\. (4.31) 

JA'B' V £4 £3 ^ J ^ JS 

Let us now assume, instead, that A is still not vertex of squares in Ai or A 2 , and 
that B is vertex of some square in Ai. We can then dehne S and C as before; this 
time, g^ in the segment BC is not dehned as the reparametrization of the image under 
/ of some segment xbXc, but as the affine interpolation satisfying g^i^B) = f{B) and 
gz{C) = f{C). However, the estimate (4.1) and the property (4.2) immediately imply 


that, exactly as before, the curve 70 intersects the image oi BC under g^ (which is the 
segment f{B)f{C)). If we call P the hrst point of intersection (starting from /(a:'^)), we 
can argue exactly as before: we dehne 7 as a slight modihcation of 7 , which is this time 
the union of the curve 70 from /(x'^) to P with the segment PB = Pf{B). In this case. 


instead of (4.30) we get the estimate 


^'(7) < f~ \Df\ + f{B)f{C) = f \Dg,\ + f{B)f{C ), 

Jx\x'^ JA'B' 


and since by the estimate of Lemma 4.3 we have of course 

2 


f{B)f{C)<- f \Df\, 
r Js 


then also in this case we get the validity of (4.31): of course even the better estimate 


without the big term in parentheses is true, but it is simpler to consider the same esti¬ 


mate (4.31) in both cases. 


Let us hnally consider the general case for the segment AB: both the points A and B 
can be vertices of some square in Ai or in A 2 . Nevertheless, as noticed above, in the cases 
already considered the path 7 coincides with 70 from the starting point x'j^ to almost the 
hnal point and only a small part near the end has been modihed. As a consequence, 
it is obvious how to deal with the case when both the points A and B are vertices of 
squares in Ai or in A 2 '- we let 7 be the path which coincides in a large central part with 
7 o, and we apply one of the above described modihcations both near the starting point. 
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and near the endpoint. Of conrse, for a side AB where we have done two modihcations, 


instead of (4.31) we will have 


n\^)< f\Dg2\ + (20A0"^ + 22"^ + ^eX+‘^f \Df\, (4.32) 

Ja'b' V ^4 Ss 5 J r J-Rab 

where as nsnal TZab is the nnion of the sqnares having either A or i? as a vertex, which 

contains both S and the corresponding sqnare aronnd A. In this way, we have hnally 

dehned in the whole grid Q, and by constrnction it is clear that the map g^ is injective 

and coincides with / on dfl. We conclnde this step by evalnating the integral of \Dg^\ on 

the bonndary of the different sqnares: for any sqnare S of the grid, we call again 5"^ the 

intersection with O of the nine sqnares aronnd S. 

If 5 G we do not need any particnlar estimate, in the next steps we will only need 
to nse that g^ coincides with ips on dS. If 5 G A 2 , instead, we know that g^ coincides 


with gi on dS, hence we only need to keep in mind the estimate (4.28) already fonnd in 
Step III, which (adding on the fonr sides of S) gives 


[ \Dg3-M-u\dn^ < ( 

JdS V 


< ( 4080 — + 44 — + -e 5 )r < k( — + — ]r 
£4 £35/ V ^4 £3, 


WSeA2, (4.33) 


Finally, if 5 G ^3 U A^, then J^^lDg^l is the snm of the integrals on the fonr sides; for 
each side A'B\ either gs = g 2 , and then of conrse 


X\Dgs\dn^= i\Dg 2 \dn\ 

JA'B' JA'B' 


(4.34) 


or \Dg^\ = for some curve 7 = 'y{A,B) dehned as above. Since the length of 


7 can be taken as close as we wish to the length of 7 , from (4.32) we derive 

\Dg2\ 


^\Dgs\< 

JA'B' JA'B' 


f 2050 — + 23 — + £ 5 V+ - 




£4 


£3 


r Js+ 


\Df\- 


(4.36) 


Since (4.34) is stronger than (4.35), we get the validity of (4.35) for any side AB of any 


square 5 G ^3 U A^. As a consequence, adding (4.35) on the four sides and keeping in 
mind (|4.29|) and the fact that £1 -C 1, we get 


!^\Dg^\dH^< (^\Dg2\dH^+{ 

JdS JdS V 


JdS 

< ^ / m 

Sir JS+ 


8200 — + 92 — + 4£5 


I dS+ndfi 


\Df\ 


£4 


£3 


r + 


12 


JS+ 


\Df\ dV? 


8200( - + -+£5 
+4 £3 


(4.36) 


V iS G A^ U Aa 


Step VI. The piecewise linear function g^ : Q ^ . 

In this step, we want to dehne a piecewise linear function g^ : G ^ this one will 
hnally be the correct map on the grid Q, in the sense that our approxim ating function 
will coincide with 5 f 4 on To do so, it is enough to apply Proposition 4.17 to the map 
gs and call g^ = g the resulting map. As a consequence, g^ is a piecewise linear map on 
Q, which coincides with / on dQ. Moreover, since for every side AB of the grid the map 
g 4 on A'B' is an interpolation of values of g^ on A'B', of course we still have g^ = ips 


on the boundary of every square 5 G Mi. Moreover, (4.17) and (4.18) imply that the 


estimates (4.33) and (4.36) are valid also with g^ in place of g^. 
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Step VII. Definition of the approximating function /^ : —)■ 

We are almost at the end of the proof, since we can hnally dehne the reqnired fnnction 
/e- We set /e = (74 on the grid Q, hence in particular /^ = / on dfl; by construction, 
is injective and piecewise linear on Q. To keep the injectivity, is then enough to extend 
/e in the interior of any square S, in such a way that fs remains continuous and injective 
on it. We will argue differently on the different squares. 

If 5 G ^ 1 , then we know that fs = (fs on dS, so we extend /g = ips on the whole 


square S: this map is continuous and injective by construction, and by (4.1) of Lemma 4.3 
we know that 


\Df,-Df\<r% VSeAi. 


(4.37) 


Let us now consider a square S G A 2 '. in this case, we want to apply Theorem 3]^ to the 


function (p = (74 on dS. Keeping in mind the generalization of Theorem 3.1 observed in 


Remark 3.2, and calling /g the obtained extension on S, we get 


[ DfJx) — M dx < Kr [ \Dg 4 ^{t) — M ■ T{t) dVfit) 

Js JdS I 


(4.38) 


as soon as the estimate 


[ Dgi{f) - M ■ r(f) d'H}{t) < r(5MAx||M|| 
JdS 


holds. Thanks to (4.33), which holds also with 774 in place of 773 as noticed in Step VI, and 


recalling that ||M|| > £3 by dehnition of A 2 and that 5 max is a small purely geometric 


constant, the latter estimate is true thanks to (4.22). As a consequence, recalling that 


5 is a Lebesgue square with matrix M and constant 5, that 5 -C £5 by (4.23| and by 


dehnition of <5, and using also that £4 < 1, from (4.38) and (4.33) we get 


j \Df,-Df\< j \Df,-M\+ j \Df-M\<Kr^{"-^ + "A G ^2 , (4.39) 

Js Js Js \£4 £ 3 / 

where K is as always a purely geometric constant. 

Finally, let 5 be a square in ^3 U ^ 4 . This time, since 774 is piecewise linear on dS 
and S is (a 2-biLipschitz copy of) a square of side 2r, we let /g on S be the extension 
of 5^4 given by Theorem 2T, keeping in mind also the generalization of Remark |2.2 The 


estimate (2.13), together with (4.36), which is valid also with g 4 ^ in place of by Step VI, 
gives then 

L\DM<Kr f^\Dg,\dn^ < - f \Df\+Kr f |Z}/| rf-Hi + ^ ^ 

Js JdS El Js+ JdS+ndn V ^4 ^3 / 


Since 


/ iD/g - Df\ < / |D/g| + L\Df\ < / |D/g| + / \Df \, 

j S 'J S *-'tS 'J S 


fs JS Js JS Js+ 

since £4 < 1 , and since iF is a purely geometric constant while £4 < 1 , we deduce 

L\Df,-Df\<^f \Df\ + Kr f + + G A • (4.40) 

Js £1 Js+ Jds+ndn \£4 £3/ 
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Notice that the same estimate holds true also for S G ^ 3 ; nevertheless, in this case we 
can say even something more. Indeed, since 5 is a Lebesgue square with matrix M = 0, 
by Dehnition ED we know 


Js+ 


\Df\ = \Df-M\< 36rH < 36r% . 


IS+ 


As a consequence, for squares in ^3 we can deduce from (4.40) 

L\Df,-Df\<Kr [ + + - + V 5 g^ 3 . (4.41) 

•jS •JdS+ndU V^4 ^3 ^1/ 

Step VIII. Conclusion. 

We are now ready to conclude. Indeed, by construction is a hnitely piecewise affine 
homeomorphism which coincides with / on dQ. Moreover, it is immediate by construction 
that ll/e - /IIloo and ||/e - /|| are as small as we wish (it is enough to have chosen at 
the beginning r small enough); as a consequence, we can assume that they are smaller 
than e/4 each. As a consequence, to get ( |4.19| ) and conclude, we only have to check that 

l|B/,--D/lli. < h (4.42) 


Thanks to (4.37), (4.39), (4.40) and (4.41), and calling, for j = 1, 2, 3, 4, 


Hj — UI iSj ! iSj G Ml 




we have 

[\Df,-Df\=[ \Df,-Df\+[ \Df,-Df\+[ \Df,-Df\+[ \Df,-Df\ 


</TA+^ + £5)|SJ|+ ^ IB/I+ ^ Kr 

\e4 £3 EiJ r.Si&A4 iicS.-pAUJA. 


i-. 


dstr\an 


\Df\. 


(4.43) 


Let US now recall that, for each square Si of the grid, the set 5/ is the union of the nine 
squares around it (to be precise, to those which belong to hi). As a consequence, calling 
for brevity At = also recalling ( |4.25[ ) we have \At\ < 9| ‘5/ < £ 2 - 

K £ie £ 


Thus, by (4.20) we can write 

K 

i-.SiGAi 


E //jB/l<9 
ItVa £^ Js+ £1 JAl 


£i 54K 6 

Analogously, each side in dfl can belong to dSt for at most eleven different indices i with 
G M 3 U M4, hence by (4.24) we get 

5: Krf \Df\ < llKr f \Df\ = llKrP{f{n)) < ^ . 

rV .,. Jds+nan Jan ' 6 


i: Si&AsUAA 


Finally, by (4.22) we have 


K 


Vir4 £3 £ 1 / 


£ 

< - 
- 6 


Inserting the last three estimates inside (4.43) we get (4.42), so the proof is concluded. □ 
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The above proposition shows that, under stronger assumptions than in Theorem 1.1 


we can obtain something better than what claimed in Theorem 1.1 Indeed, if Q is an 
r-set and / is piecewise linear on dQ, then we do not only get simply a countably piecewise 
affine approximation, but a much better finitely piecewise affine one. We can now give 
the sharpest possible result of this approximation, that is, we can prove the existence of 
a hnitely piecewise affine approximation with the weakest possible assumptions. 

Theorem 4.20. Let hi C be a polygon and let f G iy^’^(r2, M^) be a homeomorphism, 
continuous up to the boundary and such that f is piecewise linear on dfl. Then, for 
every e > 0 there exists a finitely piecewise affine homeomorphism /^ : hi —)■ such that 
Wfe - fWwn + Wfe - /||l°° < and f^ = f on dLl. 

Proof. Since hi is a polygon, there exists an r-set hi and a finitely piecewise affine home¬ 
omorphism 4) : n —)■ n. There exists then some constant H = H{Q) such that 


|D4)(x)| <H, 


det D^{x) > — , 
H 


for almost every a; G hi. Let us then define f:il^Mfasf = fo by construction, 
/ belongs to lT^’^(r2,M^), and it is cont inuou s up to dfl and piecewise linear there. As a 
consequence, we can apply Proposition 4.19 to / in hi, finding a finitely piecewise affine 


homeomorphism /g : hi —)• which coincides with / on dfl and such that 

II/.-/II, 


IVLi’i(O) 


+ ll/e - 


lL°°(f2) — PJ2 


(4.44) 


We can then define /g = /g o $: this is a finitely piecewise affine homeomorphism, of 
course it coincides with g on df, and we have 

||/e-/IU-(0) = ||/e-/|lico(H)- 

By a simple change of variable argument, we obtain then 


ll/e-ZIUpo) = / fe{x)-f{x) 


fe{y) - f{y) 


dx= fs{<^{x)) - f{^{x)) 

J 

dy < H\\fe - 


dx 


iLpO) 


detD^(^-^{y)) 

and similarly 

WDfe - ^/lUpo) = \D(h ° - D(f O <l.)(a;) 

= [Dh{Hx)) - Df{^x))) ■ Dd>{x) dx<H^Df,-Df\\ 


dx 


Li(0) • 


Inserting the last three estimates in ( 4.44[ ), we conclude that /g is the desired approxima¬ 
tion. □ 


Remark 4.21. It is immediate to observe that the assumption of Theorem f.20 are sharp. 
Indeed, assume that a homeomorphism f G admits a finitely piecewise affine 

approximation /g. Since /g is finitely piecewise affine, it is defined on a polygon, hence hi 
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must be a polygon. Similarly, is piecewise linear on dfl by definition, and since fe = f 
on dfl the same must be true for f. 


We are now in position to give the proof of our main result, Theorem 1.1 


Proof of Theorem 1.1 We will argue in a way quite similar to Proposition 4.19, we only 
need some additional care to reach the boundary of fl. First of all, we look for a piecewise 
affine approximation, the smooth one will be found at the end. 

We again start by selecting the small constants ep hrst of all, we let Si be a small 
geometric constant, say ei = 1/10. Then, since / G we can select a constant £2 

such that 


\Df\< 


eei 


MA<PfL: |Al| < £2 


(4.45) 


The next step is to write as a countable union of r„-sets. More precisely, we can take 
a sequence of constants Tn —)■ 0 and a sequence of disjoint open sets CC fl, in such a 
way that each is an r„-set, the union of the closures fin is the whole fl, and for each 
n G N we can divide the boundary of fin in two disjoint parts dftn = d~ftn being 

c?“f2i = 0, d'^ftn = d~ftn+i VugN. 

Since / G we can select these sequences in such a way that 

[ \Df\ = P{f{fn)) < +00 Vn G N, (4.46) 

and we can also take fli large enough so that that 

esi 


/o\Oi 7 /a 


(4.47) 


Naively speaking, the idea now is to try to work on each set fin separately. However, since 
/ is not necessarily piecewise linear on the boundaries of the sets we cannot simply 
rely on Proposition 4.19 for each moreover, since H has not necessarily hnite area, 
estimates like (4.37) or (4.39), where the area of a square appears, are not acceptable 
because they could give an inhnite contribution after adding. 

Let us now concentrate on Hi in order to select the constants £ 3 , £4 and indeed, we 


will use these constants only inside Hi. Arguing exactly as in the proof of Proposition 4.19 
we hrst let £3 be a constant such that 


<|. (4.48) 


ja; G Hi : 0 < \Df{x)\ < £3 or \Df{x)\ > — or 0 < det(H/(a;)) < 
then we let again 

= |£3 < ||M|| < — , det M > £3! , = |£3 < ||M|| < — , det M = o| , 


then we let £5 £4 £3 be such that 

^ ^ ^ < g 

g4 gs gi ~ 8iF|Hi| ’ 


gs <C £4^3 , 


£4 -C £ 


3 ’ 


(4.49) 


and hnally we let 6 = 6 be given by 

5 = min |^(M, £5), M G Ad"*" U Ad°} , 














54 


STANISLAV HENCL AND ALDO PRATELLI 


where 6 are the constants of Lemma 4.3 The last thing we have to hx, is the hnal value of 


the constants indeed, each is an r„-set, but then it can be regarded as a r„/iL„-set 
for every constants Hn G N. As a consequence, we can now decrease the values of 


(without changing the sets of course); in particular, also thanks to (4.46), we can 
assume that ri is so small that 


riF(fii) + 


{a: G ffi : r{x, 5) < ri} 


< 


g 2 

180 ’ 


riP{/(SJi)) < 


mK ’ 


(4.50) 


being r(a;, 5) the constants of Lemma 4.6, while any other is so small that 


rn < Tn-l , 


/ dO-n 


\Df\< 


88K ■ 2^ 


<C dist(c? Vn>2; (4.51) 


notice that the last requirement basically means that the “thickness” of any is of 
several squares. 

Having hxed the sets Qn and the corresponding r„, any is divided in a hnite union 
of squares, all with side 2r„. Let us enumerate them by saying that the squares of the 
grid of fin are iSf with I < i < N{n); then, we subdivide the squares of f2i in four groups, 
namely, A\, A], Al, Al, as follows, 

-41 = {Sj CC 111 : Si is a Lebesgue square with matrix M} G and constant (jj , 

Ai = {y CC 111 : Si is a Lebesgue square with matrix Ml G and constant (5| , 

Al = {5‘ CC : Si is a Lebesgue square with matrix Ml = 0 and constant (5| , 

.4l = {5y S‘^.4lu.45u.4j}. 

We immediately record that, exactly as in Step I of the proof of Proposition |4.19 


from (4.48) and (4.50) it follows that 


|J{5i G ^ 4 }| < riP{fli) + 4^^ {x G 111 : r{x, (5) < n} + 


< 


£2 

9 


(4.52) 


For any n > 2, instead, we simply let A!l be the collection of all the squares Sp of the 
grid of while A^, A^ and A'^ are empty. 

Let us notice that the only assumption which is true for in Proposition |4.19 and 
may fail now for the generic fin is the following one: / is assumed to be piecewise linear 


on dfl in Proposition 4.19 while / needs now not to be piecewise linear on the boundaries 
dfln'i on the other hand, this assumption has been used only in Step VI of the proof of 


Proposition |4.19[ As a consequence, we can repeat verbatim all the arguments of Steps II, 
III, IV and V of that proof for fli: so, we discovery hrst that squares in A\ and A\ can 
never touch, then we dehne a temptative modihed grid Qj with a function gl : 
then the correct modihed grid with the function gl : ^ , and hnally the function 

gl : ^ By dehnition, the function g^ is injective and coincides with / on dfli (this 
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was explicitely decided in Definition 4.15), and moreover we have the estimates 


las 


fas \^9\ ' 
.\Dgl\dH^< 
\Dgl\dH^< 


gl = (ps on dS 
M -uidn^^ k(— + — 


n 


las 


las+nani 


K 


|C/| + /(55+'^ + ^)r. 


f \Df\+ [ + + 

Js+ Jas+nani \e 4 e^J 


G 

V5 e , 

ys eAl, 


(4.53) 


BiTi Js+ -ids+nani 

where K is a. purely geometric constant (it suffices to take K = 9000 here). 

Let us now consider for any n > 2. In this case, the situation is much simpler 


than in Proposition 4.19 indeed, by definition we only have squares in A2, so we do not 


need the arguments of Steps II and III and we can directly start with the analogous of 


Step IV, which immediately gives us a function g^ 
of ( 4.29| ), that is, 

400 


Qn —t satisfying the analogous 


lDg^ldn^< 


\Df\ dW + 


IDfidn^ 


V5 G 


(4.54) 


las' ""' EiTn Js+ Jas+nan„ 

Again since there are no squares in Ai, A 2 and ^ 3 , we do not even need the argument 
of Step V, and we can simply set the function g^ = g 2 - 

Let us now observe that every function g^ coincides with /, by construction and by 


Definition 
we set g^ : 


4.15 


on d^n- As a consequence, if we call Q the union of all the grids and 


^ as 5^3 = g^ on each then the resulting function g^ is also injective. 

The last thing we have to do, before having the right of treating each separately, 
is to modify g^ so to become piecewise linear on the boundary of each we will do that 
by applying Proposition 4.18 More precisely, for every j > 2 we define O the union of 
all the segments of and Qi, 


and Cq = d 


D,-. We apply then Proposition 


4.18 


with 


C = Ci, Cq = Cq, T] Tj < rj_i, and with the function g = gz] thus, we get a function 
piecewise linear on d~Qj, which coincides with g (hence, with gz) out of a 77 -neighborhood 
of d~Qj. As a consequence, the function gj is different from gz only on the boundary of 
squares which meet d~flj. Then, let us define the function ^3 : ^ as ^3 = gj on 

the boundary of squares touching d~Qj, and gz = gz on the boundaries of all the other 
squares. By construction and by Proposition |4.18 the function gz is injective, piecewise 


linear on each dilj, and the estimates (4.53) and (4.54) are valid with gz in place of gz- 
Now, for every n G N we apply Proposition 4.17 to the set fin with the function gz on 
and we get a new function g^, piecewise linear on and coinciding with gz on dfln - 


Finally, we define g^ \ Q ISr as g^ = g2 on every Q'^: also this function satisfies (4.53) 


and (4.54), and it is piecewise linear on the boundary of each square of any grid. 


We are now ready to define the piecewise affine approximation fp. indeed, for every 
n G N, the function g^ on is injective and piecewise linear on the boundary. Exactly 
as in Step VII of the proof of Proposition 4.19, we can then define /” on fin, which is 
a finitely piecewise affine function coinciding with g^ on dfln, and then set /^ : D —)■ 
as the function which coincides with /” on every fin- This function is by construction 
a countably piecewise affine homeomorphism, and also locally finitely piecewise affine; 
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moreover, it is clear that ||/ — fe\\L°°{Qn) and ||/ — /£||li(o„) are as small as we wish as 
soon as the constants r^, have been chosen small enough; in particular, we can think that 
both them are smaller than e/4, and moreover we get the fact that /^ = / on dfl whenever 


/ is continuous up to dil. As a consequence, to conclude the proof of Theorem 1.1 for 
what concerns the piecewise affine approximation, we just have to check that 


m, - B/|U.(n) < 2 


(4.55) 


This will be obtained arguing almost exactly as in Steps VII and VIII of the proof of 


Proposition 4.19, More precisely, let us start with fli: repeating verbatim the arguments 


bringing to (4.37), (4.39), (4.40) and (4.41), this time from (4.53) we get 


[ \Df,-Df\<rls, 
Js 


I \Df, - Df\ < KrU ^ ^ 

/s V^4 ^3 


pDf,-Dfl<Kr, f + + ^ + 

Js JdS - - 1 ^ ^ „ I 


/a5+naoi 


£4 £3 £1 


G Ml, 
e Ml, 
V5 gM1, 


i\Df,-Df\<-j \Df\+KrJ \Df\dV} + Kr-^ + ^Arl ^SeA\ 
Js £1 Js+ JdS+f\dQ.i \£^4 ^3/ 


4 • 

(4.56) 

Instead, for every n > 2, we apply Theorem 2T -recalling also Remark 2.2- to the generic 
5, which is a 2-biLipschitz copy of a square of side 2r: then, also from (4.54), we get 


i\Dh-Df\< i\Df\+i\Dh\< f \Df\ + KrJ^\Dg,\ 
Js Js Js Js+ JdS 


< 


/s 

K 


f |Zl/| + Ar„/ 
Js+ Jdi 


\Df\ 


(4.57) 


V5 G M^ 


£1 Js+ Jds+ndn^ 

Notice that this estimate is better than the corresponding one for Proposition |4.19 


namely, (4.40): indeed, there we had also the additional term K{e^/ei + ei/e^)r‘^, which 
now would be quite a problem since in principle may have inhnite area. The reason 
why we do not have this term now, is that it was coming from the interaction between 
squares in M 4 touching squares in Mi or M 2 , while now in we only have squares in M 4 . 


The very same argument as in Proposition 4.19 implies again that every square can 
belong to at most for 9 different squares S, and that every side of some can 
belong to dS^ fl dVLn at most for 11 different squares of the grid of as a consequence, 
adding (4.56) and (|4.57) for all the squares of the different grid, we hnd 




\£4 ^3 


£ 

< - 


£ £ 
X + X 


+ -)vil 

K f 
79— / 

£1 JQ,\fl 

HIKY. 

n£N 


+ 9 


K 

£1 


\Df\ 


\Df\ + nKY: 


nSN 


' dUri 


\Df\ 


8SK ■ 2^ 


£ 

< - 
- 2 


where M 4 ’'" is the union of the sets for all the squares S G M^, so that by (4.52) 


we have IM 4 | < 9| U {Si G Al}\ < £ 2 , and where we have used (4.49), (4.45), (4.47) 
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and (4.51). As a consequence, we have established the validity of (4.55), so the proof of 
the existence of the required piecewise affine approximation is concluded. 

As already remarked, once the piecewise affine approximation is found, the existence 
of the required approximating diffeomorphisms is exactly the content of m Theorem A], 
so we have finished our proof. □ 
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